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SYNOPSIS 


Composite materials wilh then well known advantages have tound 
widespiead use in various enginccnng sUuctiues. These have high 
sfiength to weight lalio, stittness to weight laliO/ excellent fatigue re- 
sistance and much better theimal expansion chaiacteiistics than metals 
They outpeifoim the conventional stiiictuial mateiials in many aieas 
Multi-ply laminates employing continuous fiber composites aie the most 
impoitant class among composites for structural applications 

Hxeicise of complete control on all aspects of any manufacturing / 
labncalion pioccss is veiy difficult This leads to unccrlainlies in the 
maleiial pioperlics and geometric dimensions of the components This 
IS especially true foi laminated composites because of the laige numbei 
o( parameters that aie associated with it's fabrication. Uncei lainlics in 
several factors like the fiber volume fiaclion, fiber orientation, volume 
fiaction, fibre-matrix mtciface parameteis etc. are inheienl in fibci 
icinforced composites. The dilference between minimum and maximum 
values of some of the elastic pioperlies can sometimes be as high as 
iOO percent The uncertainties in the faclois mentioned above aic m 
turn leflectecl on the chaiacleristics of the lamina slilfness parameteis 
like the longitudinal modulus, transverse modulus, Poisson's ratio, shear 
moduius etc. These are treated as the primary vaiiables, as these are the 
basic parameters that are considered for formulating a structural analysis 
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pioblem When painmeteis like modulus, density, Poisson's lalio etc 
aie random, the deiivcd lesponse chaiacleiislics like deflections, nahiral 
liequencies, buckling loads, sliesses, stiains etc aie also landoin, being 
functions of these basic random mateiial paiameleis 

Conventional methods of analysis of stiuctures assume the vaiious 
paiameteis involved in the analysis — like system chaiacteiistics, bound- 
aiy conditions etc — to be deterministic In actual piactice, howevei, 
these involve mheient unceitainties and aie seldom dcteiniinistic. In 
most cases these paramcteis can at best be expicssed in a piobabilistic 
sense only Use of the deterministic techniques yield the solution toi 
the mean lesponse only, which may load to unconseivative lesults, For 
actospace and other sensitive applications it may be necessary to adopt 
the moie accuiate probabilistic appioach of structuial analysis 

The volume of literatuie that exist for problems related to isotropic 
metallic strucluies with unceitainties in paiameteis is considerable Lim- 
ited literature is available on response for stiucluies made of composite 
materials, involving randomness in system properties The ptesent woik 
was motivated by the fact that the literature available on Iho study of 
composites, considering the system properties as random vaiiables is 
limited and there is a need for syslematic investigation of the effects of 
this randomness. 

The general formulation /solution techniques employed in this thesis 


are; 
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1 Peilurbation 

2 Monlc Cailo stmulalion 

III pciluibation loiimilation the landom painmeteis involved aie o>c- 
picbsed as a truncated asymptotic expansion, by intioducing a siuall 
jHvnnu'fcr This kind of appioximatioii can be used when the dispci- 
sion IS binallei than the mean value of the random vaiiable This is 
true for many materials and the appioach can be applied to most of 
the pinctical engineering pioblems Ihe system equations aie analyzed 
to obtain relations in the deterministic mean and the supeiimposod 
zero-mean inndom part. These lelations aie solved to tind the moan 
and the variance of the response 

In general, the governing equations of any stiiictural analysis prob- 
lem, with the system parameters modelled as random variables, can be 
expressed as. 

wheic [L/jj is the random system matrix, is the random response 
vector and {(/,} is the forcing vector, Heie L^j aie known functions of a 
set of primary RVs 6/^ Af are unknown and also functions ot bj^ The 
character of the pumary random variables is assumed to be completely 
known. The forcing vectoi {qt,} can be random or deterministic in 
nature. The problem, therefore, is to find the statistics of when the 
statistics of the primary RVs, bj^ are known 
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Hluo the ivindomnefab in system piopeilies is taken into account 
using a peituibntion technique Any random vaiiable can be exptessed 
as a zeio mean random vaiiable supciimposed ovei a mean. This can 
be lepresented as 

randovi variable = uicau + ^cro mean 7 andorn vaiiable 

For many piaclical engineeiing situations the random component is 
small compared to the mean By mtioducmg a small per luiba lion 
paiamelei i, the above equation can again be expressed as 

= IiV'‘-\-eRV’ 

Based on those, equations to solve for the mean and standard deviation 
(SD) ol the unknowns like plate deflection, naluial frequencies, initial 
buckling load etc are developed From these the required unknown 
parameleis are computed The above technique is general in naluic 
This can be used to solve a variety of pioblems lelaled to analysis of 
structures, involving random system parameters. 

In Monte Cailo simulation, a random sample of the system pa- 
rameters is generated first. This is used to gel a sample of the desired 
results, by repeated numerical experiments With the results from many 
such samples, statistics of the unknowns can be found. In the present 
work this technique is used for validation of results obtained with the 
perturbation technique. 
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The analysis technique is applied to composite laminates, based on 
the Classical Laminate Theory (CLT) Heie the basic elastic piopeities 
ol a composite lamina — like elastic moduli, Poisson's latio etc — 
aie assumed to be independent landom variables Starting fiom these, 
expiessions for the expectation of various intermediate teims — like 
the leduced stiffness matrix Q^j, bending stiffness maliix etc — 
aie obtained The above equations, aie then utilized to solve tot the 
charactciislics of the random i espouse using the peiturbation technique 
outlined Several application piobicms, employing the pioposed method 
are discussed The basic formulation of these problems is done by eneigy 
methods or other suitable conventional techniques available A biief 
examination of the elfect of coupling of the primary random vaiiables 
is also done. 

The fust problem studied using the technique developed is deflec- 
tion of composite plates with random system propeilies subjected to 
transverse static loads The basic formulation of the problem is clcvol- 
opecl based on the CLT. The solution of the equations lesulting fiom 
the foimulalion using energy methods is attempted by Rayleigh-Rilz 
technique. Rectangular plates with specially orlholropic and mtdplane 
symmetric stacking sequences are considered. Parametric study with 
variation in stacking sequence, fiber orientation, boundary conditions, 
aspect ratio, material properties etc. are presented. The approach and 
some assumptions used are validated by Monte Carlo simulation. 
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Some sludieb lolatcd to fiee vibration of rectangulai composite plates 
with randomness in system paiameteis are attempted The effect of 
landomnesb in system paiameters is intioduced using the peiturbation 
technique The governing equations obtained using energy formulation 
aie solved by the Rayleigh-Ritz technique The validation of the results 
obtained is done by Monte Carlo simulation The statistics of landoiii 
natural liequencies are obtained for vaiious types of plates A paiametric 
study changing boundaiy conditions, aspect ratio, slacking sequence etc 
lb also presented , 

Discussions on problems in the aiea of initial buckling analysis ol 
rectangular composite plates with randomness in system parameters 
are also presented. Here again the basic formulation is based on 
the CLT The randomness in system parameters aie analyzed using 
the peiturbation technique. The goveining equations aie obtained 
using energy formulation The solution is attempted by Rayleigh-Rilz 
technique A parametric study for investigating the effect of SD of 
aspect ratio, slacking sequence etc on the initial buckling load is also 
conducted 

Lastly some of the general observations and conclusions are pre- 
sented. Most of the results indicate a linear trend for the response SD 
in relation to variation of SD of input RVs. This indicates the negligible 
effect of higher order terms in the approximating functions used For 
the problems considered the mean of the derived parameters coincide 
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wilh the values obtained using the usual deteiministic analysis Com- 
paiibon ol the piesent approximation with the results obtained fiom 
Monle Calk) sinuilalion shows that for the range of SDb of input RVs 
consideicd, accuracy of lesults fiom the piesenl formulation is excel- 
lent At low levels of SD of input RVs, change in fibei orientation, 
aspect latio etc do not have much influence on the output SDs The 
nature and magnitude of influence of each one of the primary RVs is 
dilteicnl If we consider umdiiectional laminae, variation m Eu has 
maximum influence on the SD of results The SD ot buckling loads and 
natural fiequencies are stiongly dependent on the mode shape undei 
consideiation 

The analysis biings out clearly the impoitance of considering the 
randomness in material properties The classical deterministic analysis 
may at times lead to nonconservative estimates. The amount of erioi, 
howevei, depends on the dispeision of the mateiial parameteis The 
sensitivity of the SD of the response is a complex function of many 
factors like the basic random variable unclei consideiation, aspect ratio, 
stacking sequence etc. 
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Chapter 1 


INTRODUCTION 


Reseaich on new materials is an exciting, rapidly changing field This is 
mainly because of their applications in aeiospace industries which de- 
mand belter mateiials New alloys, plastics, composites etc are always 
being presented as candidate mateiials These materials, because of 
their supeiior poiformance in terms of their sliength, stiffness, strength 
to weight ratio, stiffness to weight ratio, tailorable properties etc throw 
open new possibilities and challenges for the analyst/designei Fiber 
reinfoiced and advanced composites form an impoitant pait of this 
new challenge 

1 1 FIBER REINFORCED COMPOSITES 

Composite materials may be defined as that class of materials in which 
there are two or more chemically distinct constituents on a macro 
scale, having a distinct interface separating them A bioad spectrum of 
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materials come under this definition 

Polymer matrix composite mateiials with then well known advan- 
tages at noimal tcmpeiatures have found widespicad use in vaiioiis 
engineeiing stiuctuies In view of the fact that composites have high 
strength to weight ratio, stiffness to weight latio, excellent fatigue lesis- 
tance, and much bettei thermal expansion chaiacteiistics compared to 
metals, they outperfoim conventional stiuctuial materials 

The fundamental difleience between composites and conventional 
stiuclural mateiials is that composites aie basically heleiogcneous matc- 
iials, though we treat them as homogeneous at macro level This is an 
added advantage foi the designei in terms of optimization of strucluial 
response The design of the slructuie can be interlinked with the de- 
sign of the composites used, to produce optimized stiuctuies with much 
superior performance than structures employing conventional mateiials. 

Fiber reinforced composite materials have strucliiially superior qual- 
ity fibers bonded together by a low peiformance matiix mateiial. Re- 
inforcement material is introduced in the composite in either discon- 
linuous/chopped form or continuous form. The building block for a 
structure is a lamina. However, a lamina being very thin, cannot be 
employed for withstanding external loads. For this reason, laminae 
are stacked together to make a laminate. Multi-ply laminates employ- 
ing continuous fiber composites are the most impoitanl class among 
composites for structural applications. 
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The fibers that are used currently are glass, caibon, giaphite, aiamid, 
and boion Glass fibre composites aie widely used and their piocessing 
technology is highly developed because of low cost They are being 
used m many structuies like ship hulls, antennas, spoils goods etc 
Carbon/graplute fibers have supeiioi stiffness, better fatigue stiength, 
low thermal expansion coefficients, and good heat resistance Their 
performance is good under compressive as well as tensile loads Caibon 
fibers are being widely introduced into aerospace industry because of the 
above advantages they offci. Kevlai fibie (aiamid) composites have veiy 
high tensile strength combined with low density, high toughness, and 
high tensile modulus Howevei, they have low compiessive stiength 
In view of this, they aie not suitable in situations where compressive 
stresses aie hkely to develop. Composites with boron fibei, exhibit 
excellent compressive as well as tensile load bearing characteiislics, and 
are used primaiily for aerospace applications. In certain situations hybrid 
composites are employed to take advantage of the superior piopcilies 
of different types of fibers. In many cases, they result in weight and 
cost savings when compared to structures made up of one type of fibei 
only. 



1 2 MANUFACTURING OF COMPOSITES 


The manufactunng/piocessmg techniques employed for composites is 
veiy much dilferent from that used for conventional materials like met- 
als Normally, manufacturing of a continuous fibei reinforced composite 
structure, like a plate, involves 

• Designing the laminate layup sequence taking into consideration 
the application and the loads involved 

• Actual laying up, with the design fiber orientation and the indi- 
vidual lamina thickness 

• Curing, using a specified curing cycle, depending on the structure 
and the fiber and matiix materials involved. 

Each stage of manufacturing of composites involves processes which can 
often be difficult to control to a great accuracy lesulting m variability 
in the mateiial pioperties The main sources of this property variability 
in composites can be listed as follows [1]: 

• inherent and production-related fiber and matrix property vari- 
ability, 

• variations in intermediate materials (eg, prepregs, sheet molding 
compounds); 


variations in fabrication processes; 
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• local and oveiall vaiiations in fibei volume fiachon, 

• vaiiations in fiber oiientalion resulting fiom various souices such 
as resin flow and poor initial placement, 

• voids 

Depending on the stiucluie and the mateiials utilized, a wide va- 
riety of piocesses like hand layup, filament winding, pultrusion etc 
are employed foi the manufactuiing of composites The mechanical 
and physical piopcrties of composites depend on the basic building 
blocks — the fibei and the matrix material — as well as the various pio- 

t 

cessing related factois like fibei volume fiaction, spatial distribution and 
orientation of fibeis, void fraction, interfacial bond chaiacleristics etc In 
resin-matiix composites, the properties of matrix material aie strongly 
dependent on processing conditions [1] The influence of processing 
conditions has been brought out clearly by Maekawa et al. [2]. The 
diffeience in various properties of laminates piocessed undei similar 
conditions are shown in Figure 1.1 [2]. Fiber-matiix interface properties 
of composites are strongly dependent on the processing of fibers and 
the curing cycle of the composite. The stiength of the fiber matrix 
interface plays an important role in determining the properties of the 
manufactured composite. Thus, owing to the inherent uncertainties in- 
volved in manufacturing/processing techniques of composites, the end 
product can have significant variations in the properties about the de- 
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Figure 1.1. Influence of molding lot series on four mechanical characteristics of 
a CFRP laminate 

sign values Experimental data from various sources [1] endorse these 
observations (Table 1.1). It is interesting to note that experimental data 
from different divisions of the same company show this deviation in 
mechanical propeities. Non-standard testing methods could be one 
leason foi such property vaiiability. But even the caiefully fabncated 
materials show such property variation [1]/ indicating a more inherent 
cause. The variability in characteristics of the basic building blocks of 
laminates — the &ber and the matrix materials — may add to the dis- 
persion in the properties of the end product. There are a large number 
of parameters that control the basic material characteristics (like elastic 
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Table 1 1 , Variation of reported unidirectional properties for a widely used 
Graphite/Epoxy system, from five sources (Sources Major Airframe Company 
Reports) 
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moduli, Poisson's ratio etc.) of composites Sometimes these can cause 
considerable scatter m the obseived experimental data [1], as shown in 
Figure 1.2 The reference also gives plenty of examples of uncertainties 
in composite material propeities like tensile strength, tensile modulus 
etc. 


A set of experimental observations described by Maekawa et al. [2] 
graphically describes the kind of scatter usually observed in experimental 
data — Figure 1.3. An analysis of the above results indicate that some 
of the properties of the composite materials as a function of the elastic 
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Figuie 12 Tensile strength and tensile modulus of E-glass woven rov- 
ing/polyester hand-layup laminates 

modulus, may be closely approximaled by Gaussian distnbulions. 


1 3 SYSTEM MODELLING FOR DESIGN AND ANALYSIS 

Normally, convenlional methods of design and analysis of structures, 
models the system by assuming the various parameters involved — like 
system characteristics, boundary conditions etc — to be deterministic 
Usually during the modelling of the structure, an empirical term called 
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Figure 1.3 Correlation between tensile sire' cth and four mechanical charac- 
teristics 

the factor of safety is introduced at different stages It is expected that 
the factor of safely covers the risks associated with the actual property 
vai lability in the structure 

In practice, however, the system parameters aie seldom deterministic. 
Most of these parameters can be appropriately expressed in a probabilis- 
tic sense only Thus, we can have a different class of problem modelling 
in which the parameters involved in the analysis are either deterministic 
or probabilistic or a combination of the two. This aspect is illustrated in 
Figure 1.4 When parameters like modulus, density, Poisson's ratio etc. 
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Different Types of Problems 


System chai actenstics 

Boundaiy conditions 

j , ^Static 

Loading 

Dynamic 


Delci ministic 
\ or 

I PiobabilisLic 

Y 



Isotropic Mateiial P> E. Cl, V etc 

and/oi Geometiic 


Figure 1 4 Different types of problems, with randomness in parameters — 
isotropic plate 


aie uncertain, the derived response parameters like deflections, natural 
fiequencies, buckling loads, stiesses, strains etc. are also random, being 
functions of these basic random material parameters. Randomness in 
mateiial parameters/system parameters induce randomness m the fiee 
and forced response of structures. Depending on the sensitivity of the 
response parameter to the input random variables and the statistics of 
the input random variables, the statistics of the response parameter can 
vary. In certain structures this can cause mode localization [3]. Further, 
the risk of failure of the structures and its performance characteristics 
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can be evaluated accurately only by adopting a piobabilislic approach 
of analysis As brought out eailier, all these observations aie moie pei- 
tinent to composite structures in comparison to conventional metallic 
materials It may fuither be aigued that a simple factoi of safety teini 
introduced during the design process cannot take into account the inlii- 
cate effects of pioperty variability found in composites The allowable 
design stresses (like tensile yield stress, tensile ultimate stress etc ) of 
conventional stiuctuial materials (like aluminium, steel etc ) can usually 
be expected to be veiy close to then respective mean values This 
happens, because the scatter m these propeities is confined to a veiy 
narrow domaui about the mean value, owing to the effective control 
of the manufacturing/fabrication processes whicli has fewer parameters, 
in comparison to composites. However, in the case of composites, the 
allowable design stresses will normally be much less than the respective 
mean values, due to the widely scattered data [1]. A behavior similai 
to this can also be expected in the case of basic composite material 
characteristics like elastic modulus, Poisson's ratio etc 

Composite materials are likely to be used extensively in primary and 
secondary structures in aeronautical and space projects, like advanced 
aircraft, helicopter, space stations etc Computer simulation of some of 
the proposed configurations of such aeronautical and large outer-space 
installations often show closely packed/overlapping natural frequencies 
of some of the components. In such cases, even the slightest shift in 
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characteristics of the components can have pronounced effects on the 
response of the structure 

Many results are available on the distribution of ultimate tensile 
stiength of fiber reinfoiced composites [4, 5 , 6] Similai results on 
the variability of the elastic modulus, Poisson's ratio etc is limited 
Maekawa et al [2] gave some experimental results on the distribution 
of vaiious composite laminate properties like Poisson's ratio, strength, 
modulus etc The detailed studies they conducted bung out the effect 
of cuiing, piocessing etc on the material paiameleis of composites 
and their reliability As discussed before, uncertainties in several factois 
like the fibei volume fi action, fiber oiientation, void fraction, fibie- 
matiix interface parameters etc can have significant effect on the 
response of fiber reinforced composites. The uncertainties in the factors 
mentioned above are in turn leflected on the characteristics of the 
lamina stiffness parameters like the longitudinal modulus, tiansveise 
modulus, Poisson's ratio, shear modulus etc. These can be treated as 
primary variables, since these aie the basic parameters that are usually 
considered for formulating any structural analysis problem Different 
factors and uncertainties which can affect on the response of composites 
are illustrated in Figure 1.5. These considerations indicate the need for 
more accurate probabilistic approach in the analysis of these sensitive 
composite structures 

Apart from randomness in material parameters, there could be ran- 
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COMPOSITE LAMINATES 




Randomness in System Parameters 

1) Fiber volume fraction Vf 

2) Fiber orientation 

3) Lamina thickness t^ 

4 ) Basic fiber and matrix properties 

Ef A , V, & ,G„ , 


Figure 1.5: Various sources of uncertainties — composite plate 
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domness in loading, geometry and boundary conditions as indicated in 
Figure 1 4 All these can affect the response of the structuie 

1 4 ANALYSIS OF COMPOSITES 

Literature on the analysis of composite structures, with a probabilistic 
approach using parameter uncertamties is limited The general forinu- 
lation/solution techniques available are 

1 Perturbation techniques 

2 Monte Carlo simulation 

3. Stochastic finite element method 

In perluibalion formulation, the random parameters involved are 
expressed in the form of a truncated asymptotic expansion, by introduc- 
ing a small parameter This kind of approximation is valid only when 
the randomness is much smallei than the mean value of the random 
variable This condition, however, is satisfied in most of the engineering 
problems and hence, the approach can generally be adopted for almost 
all practical engineering situations The expanded forms of the terms 
are introduced into the system equations. The deterministic part and 
the random part of the solution can be separated fiom these expanded 
system equations. These two solutions can be combined to get the 
complete solution of the problem as the characteristics of the response. 
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In the case of Monte Cailo simulation, a random sample of the system 
paiameters is generated first This is used to calculate the response, by 
repeated numeiical experiments From the results of a laige numbei 
of such samples, statistics of the unknowns can be found The sample 
size IS dependent on the nature of random vaiiables and the problem 
considered Complex problems which are otherwise very difficult or 
impossible to solve using othei techniques can be solved using Monte 
Carlo simulation This is a computation intensive technique Usually 
the only limiting factoi foi the pioblems is the computing lesouices 
available It is noimally used as a tool foi validating results obtained 
thiough other methods 

Usually real life structures are so complex that it is impossible to ana- 
lyze them with conventional analytical techniques. More so, in the case 
of structures with randomness in their system parameters In stochastic 
FEM, the randomness in parameters is taken into account at the element 
level Itself [7] Normally, the basic foi mula lion techniques employed to 
take into account the randomness in finite element equations are the 
perturbation technique or Monte Carlo simulation or a combination of 
these two. The main advantage with this approach is that it is possible 
to handle more realistic structures 
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1 5 MOTIVATION 

The discussion in the pievious sections clearly demonstrates the impoi- 
tance of stiuchual analysis pioblems with paiameter randomness This 
type of analysis will be of much importance in the case of composites, 
where theie is a wide dispeision in structural system parameters Not 
much work has been done in the case of composites, in this aiea Hence 
the motivation foi the present study 

1 6 LAYOUT OF THE THESIS 

The thesis compuses of eight chapters The first chapter gives an 
introduction to the problem. A brief overview of the area of study, the 
issues involved and the motivation for the present woik is desciibed A 
brief review of the pertinent liteiature available is presented in chaptei 
two A description of the general methodology adopted is given in 
chapter three Chapters four, five and six discuss specific pioblems 
related to static analysis, fiec vibiation, and buckling of composite 
plates, along with the results Chapter seven describes the general 
conclusions, limitations of the present work, and suggestions foi furthei 
research. Chapter eight briefly describes a method to take into account 
the coupling of basic material parameters 
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17 SUMMARY 

The system paiameters of stuictures composed of composite mateuals 
aie often much moie unceitain in compaiison to those made of conven- 
tional materials Judging from the available test lesults it appeals that 
the leal life values have a wide dispeision Commonly, toi conventional 
methods of analysis of composites, all the parameters that describe the 
system aie considered to be deterministic Thus the results, so obtained 
do not leflect the effects of the real life unccitainty of the parameters 
involved and may pi edict the system behavioi inaccuiately Foi more 
sensitive applications like aerospace stiucluies, the problem of analysis 
of composites with uncertain paiameteis call for a probabilistic approach 
The pertinent literature available in this aiea and related fields is 


reviewed in the next chapter 



Chapter 2 


LITERATURE SURVEY 

2 1 INTRODUCTION 

A considerable volume of Iiteiature exists foi pioblems lelated to struc- 
tuies with unceilainties in geometiic and loading parameleis Howevei, 
very little literature is available for composite mateiials in this area. This 
may be due to the fact that there aie a large numbei of parameters that 
affect the response of composite Furthermoie, the interactive effects ot 
various paiameters are not cleaily understood. The available literatuie 
may be bioadly classified into two groups 

1 Problems involving structural or material unceitaintics for struc- 
tuies made of isotiopic materials 

2 Problems involving structural oi material uncertainties for struc- 
tures made of composite materials (generally anisotropic) 
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2 2 ISOTROPIC MATERIALS 

Seveial studies [8] exist for stiuctures made of conventional materials 
In these, the basic matenal paiameteis and/oi othei paiameteis of 
the stuicture under consideration (like geometiy, boundaiy conditions, 
loading etc) aie modelled to be landom vaiiables (RV) 

A general treatment of probability theory and it's applications to 
theory of stiuctures can be found in [9] This piovides detailed in- 
formation about problems related to static and dynamic behavior of 
slructuies with system paiameler randomness 

Ibrahim [3] has presented a comprehensive review of the various 
appioaches and solution techniques for stiuctural dynamics problems 
with paiametei uncertainties Many problems related to landom eigen- 
values, random responses etc. of discrete and continuous systems have 
been discussed The review also includes results of experimental in- 
vestigations and discussions on the effect of these types of problems 
on aieas such as design optimization, leliability analysis, sensitivity of 
sliuctuial performance to parameter uncertainties etc 

Vaicaitis [10] has discussed the fiee vibration analysis of beams with 
flexural rigidity (El) and mass (pA) as RVs, along the axis of the beam. 
Various types of distributions for the random parameters have been 

considered The formulation of the problem uses a two variable pei- 
« 

turbation expansion for the deflection. Natural frequency and normal 
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mode results foi a clamped beam, obtained by using Monte Carlo sim- 
ulation have also been presented The accuiacy of the method used 
has been found to increase with higher frequencies It has been found 
that for a beam with random chaiacteristics, the frequencies and noi- 
mal modes deviate significantly fiom that for a beam with uniform 
chaiacteiistics 

Astill et al [11] have studied stiess wave propagation through a 
sliucture with random geometric and mateiial properties undei impact 
loading Different conciete stiuctuies have been used as examples The 
Monte Cailo appioach has been employed to solve the problem The 
technique proposed has been demonstrated to be compatible with Finite 
Element Method. 

A scheme foi finding piobability distribution function of eigenvalues 
of simple stiuctuial systems (like bars, strings, shafts etc) governed 
by the stochastic wave equation has been suggested by Manohar and 
lyengai [12]. The stochastic averagmg method employed can take into 
account randomness in mass, stiffness and boundary conditions. 

A beam-column, suppoited by torsional springs at ends the has been 
studied by Shinozuka and Astill [13], The free vibration with random 
variations in spring constants, axial force, distributions of material and 
geometric properties have been considered. Statistical properties of 
eigenvalues have been found using a perturbation formulation These 
results were compared with results from a Monte Carlo simulation. The 
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pertuibation approximation has been shown to be acceptable foi cases 
where the statistical variation of properties have a limited lange 

A general approach for solving the statistics of the eigenvalues has 
been piesented by Collins and Thomson [14] The method can take into 
account correlation between matiix elements and can provide statistics 
of eigenvalues and eigenvectors Lineai equations foi dynamic sys- 
tems have been foimulated for finding the statistics of eigenvalues and 
eigenvectois The appioach has been verified by an independent Monte 
Cailo simulation 

The influence of damping uncertainty on the fiequency lesponse of a 
imvUi degtee-of-heedom dynamic system has been studied by Cavavani 
and Thomson [15] A statistical linear model that allows a piobabilistic 
description of the response in terms of second-oider statistics has been 
piesented A hypothesis regarding the statistical correlation of damping 
coefficients has also been presented. Validation of the piesented linear 
model, has been done by Monte Carlo simulation The technique is 
applicable to different forms of damping and other system parameters, 

Chen and Soioka [16] have studied the response of a multi degiee- 
of-freedom dynamic system with statistical properties, to deterministic 
excitation. The formulation employs a perturbation technique A nu- 
merical example has been presented to illusti'ate the determination of 
response statistics using the procedure given. To show the effect of 
system property variations, the results have been compared with the 
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response ot the system without considering stiuctuial piopeity vari- 
ations It was found that the piobabiUstic response always provides 
largei response values than the deterministic response In yet anothei 
paper [17] they have looked into the case of impulse response of an 
analogous system employing similar technique Statistical moments of 
the impulse response function have been obtained 

The free vibrations of a slender bar with slightly different charac- 
teiistics fiom those of a uniform bai has been studied by Mok and 
Murray [18] A simple, appioximate solution obtained by ignoring the 
higher order terms and their products, of the varying paiameteis like 
area of cross section, elastic modulus, density etc has been presented. 
Experimental results, piesented for fiee-free flexural vibrations have 
been shown to agree well with the theoretical predictions 

Grigoriu [19] has described methods for calculating probabilistic char- 
acteristics of the eigenvalues of stochastic symmetiic matrices The 
formulation has been obtained by developing approximate expressions 
for eigenvalues and eigenvectors, using a perturbation technique Sev- 
eral examples from dynamics and elasticity have been solved using the 
method and the results have been compared with exact solutions 
Weighted integral method has been used by Deodatis and Shinozuka 
[20, 21] for calculating the response variability and leliability of stochastic 
frame structures. Energy approach along with FEM has been employed 
for solving the problems. Illustrative examples for finding the stochastic 
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displacement field of a beam element with random maleiial piopeities 
have been given 

Buchei and Biennei [22] have desciibed a method toi the analysis of 
structural systems with landom properties subjected to dynamic loading 
They have used the weighted integral technique Mass and stiffness 
terms have been expanded in terms of weighted integrals Monte Carlo 
simulation and response surface method have been used to determine 
the thieshold exceedence probabilities of the response 

A procedure for analyzing stiuctuies with random stress-strain be- 
havioi has been presented by Millwater et al [23] The unceitainties 
m stress-strain behavior of a stiuctiue has been simulated by letting 
parameters like Elastic modulus, stress at yield etc. to be random vari- 
ables The approach has been used to evaluate structural reliability of 
complex stochastic shuctuies by employing a finite element scheme 
Chen and Zhang [24] have described a method to determine the re- 
sponse sensitivity for complex stochastic structures subjected to arbitrary 
deteiministic excitation The probabilistic problem has been converted 
into a deterministic one using a perturbation technique. The response 
sensitivity has been calculated with respect to parameters like beam 
cross-sectional area, plate thickness etc. 

A technique for computing the instantaneous transient response 
statistics of an undamped linear multi degree of freedom system sub- 
jected to arbitrary deterministic excitation, when random uncertainties 
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are introduced into the stiffness matrix, has been described by Piaslhofei 
and Beadle [25] Peituibation technique has been used to model the un- 
cei tain ties The lesponse unceitamty has been presented as a function 
ot the model uncertainty for an impulsive excitation function 

Bliven and Soong [26] have presented a method foi analysis of nahiial 
frequencies of elastic beams having randomly varying characterishcs 
The method is based on a lumped paramelei model The randomness 
may be due to mateiial inhomogeneities oi geometncal imperfections 
Statistics of nahiial frequencies of a simply suppoited beam have been 
calculated with randomly varying stiffness with position along the beam 
Zhang and Chen [27] have studied multi degiee of freedom systems 
With uncertainties in mass and stiffness matrix elements The formu- 
lation foi free vibration analysis has been carried out by using a fust 
order peituibation expansion of the uncertain paiameters Statistical 
properties of eigenvalues and eigenvectors have been found for beam 
and truss problems The technique employed is simple and general, 
and can be applied to some problems related to composite materials. 

Beams with arbitrary boundary conditions have been studied by Low 
[28] The formulation takes into account the effect from neighboring 
structures and boundary conditions, and incorporates them into the 
frequency equation Energy approach has been used. Exact roots of 
the equations have been found by a simple search method. Some 
experimental results are provided to validate the approach used. 
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Vanmarcke and Grigoriu [29] have piesented lesults on their study ol 
beams with random rigidity variations along the axis Static analysis of a 
cantilever beam with a unifoimly distiibuted load and concenhated load 
on the tip has been peifoimed using FEM Random elastic characteristics 
of each element has been represented by the local spatial aveiage over 
the element, with correlation function repiesented by a parameter called 
scale of fluctuation. Second Older statistics of deflection have been 
evaluated using the pioceduie 

A general pioceduie for formulating problems involving landom 
fields, using probabilistic FEM has been described by Liu el al [30]. The 
random field has been discretized into a mean vector and a co-variance 
matrix To reduce computational effort, the above corielated vector 
has been transfoimed into an uncorrelated vector using an eigenvalue 
orlhogonalisation The method has been applied to several problems, 
including a two-dimensional plane-stress beam bending problem An 
independent Monte Cailo simulation has been employed to check the 
results In another paper Liu et al [31] have described a finite element 
method applicable to truss structures for the determination of the prob- 
abilistic distribution of the dynamic response of truss structures owing 
to material and geometric effects 

A solution for free vibration of thin rectangular plates resting on 
non-uniform lateral elastic edge supports has been obtained by Gorman 
[32]. A method of superposition has been used. The stiffness of elastic 
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suppoits may have any desiied distiibution along the edges including 
discontinuities and local concentiations Results have been piesenled 
for squaie and lectangulai plates 

All the investigations mentioned above aie confined to elastic isoliopic 
materials 

2 3 COMPOSITE MATERIALS 

Limited lileiatuie is available for the response stiuctuies made of com- 
posite mateiials, with randomness m system piopeities 

Two papers by Mai tin and Leissa [33, 34] heat plane elasticity, static 
and dynamic analysis of composite lamina with variable fiber spacing. 
The fiber volume fraction has been assumed to vaiy according to some 
given functions. This results in variation of elements of stiffness matiix 
across the width of the lamina The plane stress problem of such 
a lamina is governed by a partial differential equation with variable 
coefficients The Ritz method has been used foi general formulation 
Seveial fiee vibration problems have also been solved The non-unifoim 
distribution of fibers is found to increase the buckling load by as much 
as 38 percent and the fundamental frequency by 21 percent. 

Nakagiri et al. [35] have studied fiber reinforced composite plates 
undergoing free vibration, with uncertainties in stacking sequence. Fiber 
orientation and thickness of the individual lamina have been consid- 
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ered as the basic uncertain parameteis Mean centered second oider 
peituibation expansion of stiffness matrix has been used for formulation 
of the pioblem The eigenvalue statistics has been obtained by succes- 
sively solving the series of equations obtained from the peituibation 
formulation The problems have been solved using FEM 

Engelstad and Reddy [36] have desciibed a piobabilistic miciome- 
chanics -based non-linear analysis proceduie to piedict the vai lability 
in the properties of metal-matiix composites. Monte Cailo simulation, 
with assumed piobabilistic distributions of propei ties— like fibei, matiix 
and mteiphase piopeilies, volume and void ratios, strengths, hbei mis- 
alignment etc — have been used to predict the statistics of resultant 
lamina piopeities 

The pioblem of maximization of critical buckling load of an angle- 
ply laminated plate, when the properties of the plate are known to 
be scattered around some values, has been discussed by Adah et al 
[37] Convex analysis approach has been employed to determine the 
material piopeities pioducing the minimum buckling load The results 
have been given for uniaxial and biaxial loadings, with various aspect 


ratios 
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2 4 OBSERVATIONS 

The peituibation method seems to be common foi analysis of stiuc- 
tures with paiameter unceilainties In general Monte Carlo simulation 
seems to be the choice for checking the validity of the proposed analy- 
sis/method. Some prefei techniques like the weighted integral method 
In some cases the above approaches aie implemented using finite ele- 
ment methods In most of the cases the solutions have been attempted 
using fist order approximations, due to the complexity of the lesultmg 
equations Introduction of randomness in system parameteis was found 
to result in significant deviations of the solutions from the mean values 

2 5 OBJECTIVES 

From the review of pertinent liteiature, it is observed that not much in- 
foi matron is available on analysis of composite beams, plates and shells, 
with parameter unceitainties. Techniques for finding the response of 
different types of composite laminates with randomness in system piop- 
eities are not available in the lileratuie In this thesis an attempt is made 
to study the response of ortholropic and midplane symmetric laminates 
under static, dynamic and inplane compressive loads In Chapter 3 
a general method that can be used for analyzing composite materials 
with parameter uncertainties, based on the perturbation technique is 


derived. 



Chapter 3 


GENERAL FORMULATION 

3 1 INTRODUCTION 

This chapter gives an account of the broad character of the problem and 
a general approach and formulation for problems in structural analysis 
involving random material paiameteis It has been discussed in previous 
chapters that it is very difficult to closely control the geometric/mateiial 
characteristics of composites For the present formulation, the basic 
lamina characteristics like the longitudinal modulus, tiansversc modulus, 
Poisson's ratio, shear modulus etc are considered to be random in 
natuie All other system parameters are heated as deterministic 

3 2 GENERAL NATURE OF THE PROBLEM 

The mathematical statement of a structural analysis pioblem usually 
involves three sets of parameters: 
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1 Terms describing the known system paiameters 

2 Terms describing the known forcing function 

3 Terms describing the unknown response function 

In general, the governing equations of any structural analysis prob- 
lem can be expressed as 

(i.;){A} = (31)- 

Where [Ty] is the lineai/nonlinear system matrix involving the known 
system parameters, {^4,} is the unknown response vectoi, and {g,} is 
the given forcing vector on the structure 

If the system parameters are random vaiiables, the above equation 
is transformed to 

[A'SKA"} = {5.} (3 2) 

Where [^Lyj is the random system matrix, } is Hie random response 
vectoi, and {g,} is the forcing vectoi iy can be expressed as known 
functions of a set of primary RVs bj^ Tlf are unknown and also 
functions of bp and, hence are random in nature The characteristics 
of the primary random variables are assumed to be completely known 
The forcing vectoi {g^} can be either random or deterministic in nature 
For the purposes of the study, it has been assumed to be deterministic 
The problem, therefore, is to find the statistics of Af when the statistics 
of the primary RVs bp are known. 
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3 3 SYSTEM PROPERTY RANDOMNESS 

The randomness in system properties are taken into account using a 
perturbation technique It is assumed here that all the primary vaiiables 
are independent random variables In achial practice this need not be 
true For example, the fiber volume fraction affects all piopeities 
However this is done to keep the resulting equations simple/solvable, 
since a system involving a number of dependent random vaiiables will 
give use to a highly coupled system of equations This is justified as 
a fiist approximation m the case of composites where the governing 
equations are much more complicated than those for a system involving 
isotropic materials 

Any random vaiiable can be expressed as a zeio mean random 
variable superimposed over a mean This can be expressed as 

random variable = mean -(- zero wean random variable (3 3) 

le RV^ = RV^ + RV' (3 4) 

The random variables m Equation 3.2, can be represented in the form 
given above Similarly the basic random variables can also be expressed 
in this manner. Thus we have 

A!^ = Af + A', 

(3 5) 


6|' = bf + bl 



32 


For many practical engineeiing situations the random component (i^V) 
IS generally small compared to the mean Equation 3 5 therefoie, 

can again be expressed as 

Lj = 

ylf = A'l + eA[ 

= bf + eb] (3 6) 


Here e is a scaling paramelei, small in magnitude No generality is lost 
m using the same value of e in all the cases The remaining facloi in 
the random term can absorb varying levels of this term Using Taylot 
senes expansion, some of the terms of Equation 3 6, can be expanded 
as 






T 




(3 7) 


rdb? 

Assuming b\ to be small in magnitude, the second and higher oidei 
terms aie neglected here This assumption sets a limit to the range of 
applications of the approach being outlined This assumption is justified 
if the randomness is small compared to the mean value Substituting 
for the terms of Equation 3.2 from Equation 3 6 we get. 


[Li, + eL:,] [a'! + eA',} = {q,} 


(3.8) 


Equation 3 8, can again be further expanded as: 


ii] {a"} +£[ l.J]{.4;} + j[i:,] {.4?} + £“[£;,]{>!(} = fe} (3.9) 
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Collecting coefficients of different powers of e fiom Equation 3 9, we 


= te) 


(3.10) 


= { 0 } 


(311) 


Confining to fiist older pertiiibation, the smaller teims of O(e^) and 
higher older are neglected here Equation 3 10^ lepiesents the deter- 
ministic part and Equation 3.11, lepresents the landom pait of the 
system of equations The Tayloi seiies approximations for and 

Llj fiom Equation 3.7, are now substituted in Equation 3 11 This 
results m 

+ = w pi2) 

Since the variables bp aie assumed to be independent of each other, 
their coefficients in Equation 3 12, can be equated for each I Hence we 




(313) 


j can be found by solving Equation 3 10 The random part of system 
of Equations 3.13, can now be solved for the only unknown 
each I Thus an approximation of the random response can be found 


Ap = Af + eA\ 


(314) 
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Heie A\ is given by Equation 3 7 Statistics of the lesponse can be 
found using Equation 3 7 and Equation 3 14 Thus the variance and 
hence the standaid deviation (SD) of the i espouse can be evaluated as 


Var(yl") - 

i 


' 3Af' 


(3 IS) 


The above technique is geneial m natuie This can be used to solve a 
variety of pioblems related to analysis of stiuctures/ involving random 
system paiameteis 


3 4 ANALYSIS OF LAMINATED COMPOSITE PLATES 

The analysis presented here is based on the Classical Laminate Theoiy 
(CLT), with the following basic assumptions 

• Plate deflections are small compared to its thickness 

• Thickness of the plate is small compared to it's other dimensions 

• The mid plane of the plate lemains neulial and is not stretched 

• Points of the plate lying initially on a normal-to-the-midplane, 
remains like that 

• Normal stresses in the transverse direction can be neglected. 

For the present analysis the elastic properties of a composite lamina 
are selected as the basic variables and are assumed to be independent 
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random vaiiables These random vanables — elastic moduli, Poisson's 
ratio etc — can be represented as 

0i - ii.,1 

bf = 

bf = Gf, 

= Gf, (316) 


If Qtj aie the reduced stiffness malnx terms, the simplified lelalion- 
ships between the reduced stiffness matrix terms and the basic lamina 
propel ties aie given by [38] 


Qn 

Q22 

Qi2 

Q 21 

Qee 

Q 55 


Eu 

(1 ~ M12/-^2i) 

E'i2 

(1 - /ii2/l'il) 

^ 21-^11 _ Ml 2 - 1^22 

(1 - ll\2}i-'l\) (1 “ M1'2/J2i) 

Qi2 

Gn 

G'ii 

Gi3 = G\2 


(3 17) 


Equations 3 17, can be individually expanded to represent them in an 
infinite series form Thus we get: 
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<5 515 = Uq — Uf cos 29 (3 20) 

wheie the invaiiant teims aie defined as lineai combinations of the 
reduced stiflness matrix leims Q^j as 

U \ = — [3Q11 + 3Q22 + 2Q22 ~ 4 ( 366 ] 

U 2 = 1[Qu~Q22] 

U-i ~ g [<3 11 + Q22 — 2(3 12 “ 4<3 cg] 

“ g(<3u + Q 22 + fiQi2 ~ 4(3co] 

^5 = o[Qll 4' Q'ii — 2(3 12 + 4(36c] 

o 

Ur - i[(3^w-Q55] (3 21) 

The Equations 3 20 can be written in a convenient and compact form 

like 

(Qvk = Elj<Um),U, ( 3 - 22 ) 

where aie known functions of Of,, the fiber orientation and aie 

the invariant functions defined in Equations 3 21. Again these Invariant 
functions can be expressed in the following form 

Uq — (3 23) 

Here are known constants 
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Z' 



Figuie 31 Nomenclature for dimensions and stacking sequence of the com- 
posite plate 

The elements ol the exlensional shttness matrix Atj, coupling stillness 
matiix B,j and bending stitiness matiix D,j [38] aie given by 

- M-,) 

O.; = iEf!,’ - 'El) *,J=l,2«nidO (3 24) 

wheie ill IS the thickness of the lamina, in the laminate given in 
Figiue 31 Foi the sake ot claiily only one inplane load liansveise load 
IS shown in the hgme ITowevei in gencial all inplane loads and^oi 
bending loads may be acting on the plate 

The solution .V, ol a strucluial analysis problem involving fibei 
leintoiced composites will in geneial, be a function ot the stittness 
matrix teims ol A,y, B,j and I>y as given by. 

.Y = (3.25) 

So, the expectation of .Y will be a function of expectations ol the stiffness 
matiix elements - ■E'lA;} and i7[Aj) This may be expressed 
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as 

E[.V| = r{B[A„],E[B„],E[D.,\) (3 26) 

Heie as a simplification^ we considei only those pioblems which do not 
involve the cross teims between diffeient stiffness matrix elements The 
unknown expectations in Equation 3 26 can be found with the help of 
the special compact foims of the equations alieady defined Thus the 
expectation of the stiifness matiix elements — E[A^j\, and 

aie given by 

m,] = 

(3-27) 

In the above equahons, the expectation of (Qy)^ is given by 

= E^,(C.„),EK1 (3 28) 

Again, heie the E\jJ^ terms - the expectation of invariant functions - 
aie given by. 

®l«.l = E,“.,E,“.,“™-E|<3u] (3 29) 

E[Qtj] are given by 

m.,] = E“o/v{n-®[(‘’i)'""‘’'’')} (3'30) 

The set of primary random variables foi the present study are the 
basic material properties — elastic moduli, Poisson's ratio etc. These are 
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defined by Equation 3 16 The deiivatives of the stiffness matrix terms 
with respect to aie required, foi solving the Equations 3 13 These 
can be defined in a form, similar to Equations 3 27 to 3 30 Thus we 
have 

dE\A,j\ 


abj^ 

dE[B,j] 

d/if 

dE[D,j] 


_ _ . / - I 2 

— }.iut '^i.-i) t — li2, 


dif 

~ ^ J ,, 2 , 2 ^ 

~ 2^k=\ rMJi v'<- '‘r-i; 








dbf^ 3^^=‘ dbl^ 

As an example foi the case of In = En we get 

UeT “ 

Dm,] __ u-^«,.. »g[(Q.Al .,., ,, . 

0El\ 2^^=' 0El\ ^ 

om„\ _ ^ 

0E[\ 3^^=' dE[\ 

In the above equations the teims — ^ aie given by 


0E[iQ„\ 


obr 




aE[u,] 


Heie the teims aie given by 

E <) 


Obl^ 


0E[U,i\ ^ Y-\() •<^G dE\(^ij\ 


and — ,^n ^ are given by 


dbl^ 


,5 


(3 31) 


(3 32) 


(3 33) 


= E,r.„/.>{nB[(kO'"'"‘’‘’']} (3 34) 

With the help of the above equations, the mean and the vaiiniice of 
the 1 espouse can be found out 


3 5 SUMMARY 

A general approach for the analysis of composite laminated 
plates with random system paiamcleis is presented This method 



can be applied to a wide vaiicty of pioblems, in static and dynamic 
analysis of composite stiiictures 

The appioach piesenled heie is used to solve some pioblems, related 
to analysis of fiber leinforced composite plates, in the next thiee chapteis 
The appioach and assumptions used aie validated by Monte Cailo 


simulations in each case 



Chapter 4 


COMPOSITE PLATES — STATIC ANALYSIS 

4 1 INTRODUCTION 

In this chapter, we consider the pioblem of a laminated rectangulai 
plate supported along the edges under transverse static loading. The 
geometry of the plate and the co-ordinate system used aie defined in 
Figuie 31 As stated earlier, we shall use the classical laminated plate 
theory foi deriving the governing diffeiential equations 

4 2 FORMULATION 

Considei a leclangular laminated composite plate subjected to a trans- 
verse load p{x,y). The governing equations of such a system in teims 
of displacements {vP, and w) and the stiffness matrix terms (A,j, 
and Dij) are given by [40] 
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d'^u^ d~ii^ d'^v^ 

, , , , d^V<^ , d^v^ „ d^w 

+ (yli2 + -^66 ) q ' Io .. + ■^•le-^rr ~ “^11“ 


dxdy 


dy-‘ 


di^ 


(B \'^B B 

~ '^■“16 c » •)■(:,.: “ + ^^66 j o rt ? ~ ^ 2 ()- 


d't^dy 


dxOy^ 


d<f 


= 0 


(4 1) 


, d^u° . , , 


„ , d^v^ , SV 
c^u; 

'^-^20 r^:.r»: .) ~ -022T: 


T, /r, n., S 

-®1C'5~T “ (-^12 + 2Z?66)‘ 


5'b’’ 


dx^dy 


dxdy^ 


dy^ 


= 0 


(4.2) 


^ <9'*i 

+ - “v-'- ■ “---^a'c^a/ ' ‘^"“a'cai/ 

„ a'^iu „ a'‘u° , a^u° 

+-022Tr-T “ -®u“?r:?r ~ '^^io^rr^ ~ + 25(3g)“ 


‘^> S .T^ 

H- 2(Di 2 + ^Ae) Q,p2ay2 + 4D26- 


- B. 


a(/‘‘ 

a^ii® 

5y^ 


2f> 


a®3 ‘“a'c^ay 
a^v® a'’u° 


dxdyr 


dx^dy 


av ^ a'’w” ..s 

,j -^22 Q .«i p(^i y) 


dxdy'^ 


dy^ 


(4 3) 


Equations 4 1, 4 2 and 4 3 can be simplified based on stacking sequence 


4 3 SPECIALLY ORTHOTROPIC LAMINATED PLATE 

In specially ortholropic laminates all the elements of [B] matilx aie 
identically zero, = ^ 2 g = 0 and Bie == B> 2 e = 0. The only nonzero 
stiffness matrix terms piesent are Bn, B 22 , Bi 2 and Bge. In view of 
this. Equations 4 1 and 4,2 are not coupled to equation 4 3. Since we 
are interested in only bending deflection. Equations for u® and u® are 


not considered 
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Foi specially oithotropic plates subjected to tiansveise loads, Equa- 
tion 4 3 reduces to the following foim 


7-. n/T-. S T-7 / \ 


(4 4) 


4 3 1 Simply Supported Plate 

The boundary conditions for deflection w, when all the four sides of the 
plate aie simply suppoited aie given by the following set of equations 


along x = 0 and x = a for all y 

w 

M, 


= 0 

d'^w d'w 


= 0 


(4 5) 


along y = 0 and y = b for all x 


V) = 0 


- , _ d'^w d'^w 

= 0 


(4 6) 


Navier's method can be applied to solve the Equation 4 4, for the above 
set of boundary conditions. A senes approximation for the deflection 
of the plate, w(x,y) is written as* 


(4 7) 
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where X„i(x) and are functions that satisty the set of boundaiy 

conditions given by Equations 4 5 and 4 6 In teims ot beam character- 
istic functions, Equation 4 7 can be lewritten as 


w(x,ij) = 


E oo ^ — ^oo 

,„=l ^„=1 


m'TT'v nTTj/ 

sm — — 

a 0 


(4 8) 


The tiansverse load distiibution p(x,i/) can also be similarly expiessed 


as 


. ^ — nDQ fivfi V 

PKv) = L,„=iLn=l‘?"*» /, 


ruTn; titti/ 
a 


(4 9) 


When Equations 4 8 and 4 9 aie substituted in the governing equation 
4 4, we get. 


Y-^oo v-'oo _ /ruTTicV „„ . /rn7rx\^ /n7rij\‘ 

) +2(D.H-2D„)(— ) (Y) 

/Tnr!/\*1 


H~ ^2'2 


Qmn (4.10) 


Equation 4 10 is of the form* 


E OO ■^—>00 _ 

,H«=1 2_yu=l 


(4.11) 


where 


Lmn = -^ll'^ni + 2(Z)i2 + 2Dca)Xl^X^^ -\- 

with X^^'^ and A» = 2f 
The expectation of L,nn is 


ir[i„J = J5[D„]Al,+2(B[D,2]Al,+2i!;Poo])A»A“ + jrp22jAl(4.12) 
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Where E[D^j\, the expectation of £>,; is defined by Equation 3 27 The 

babic landom vaiiables are assumed to be independent Here they 

are defined by Equation 3 16 Fiom Equation 4 12, taking deiivalives 

with respect to each of the basic random vauables, we get* 

dE[L„„.] _ dE[Dn]., . JdE[Di2\ . ^dE[D,,]\ ,, ,3 , dE[D,2],, 
dbl^ ~ dbi^ ^ [ dbf^ " 

(413) 

The derivatives that appear in Equation 4 13 are defined by Equations 3.31 
Thus, with the help of the expressions defined above and using 
Equations 3 10, 3 13, 3 14 and 3 15, the mean and variance of the 
response w^^{x,r/) can be found out, provided the statistics of the basic 
random variables fcf aie known. 

4 3 2 Clamped Plate 

For a plate clamped on all edges, the differential equation cannot be 
solved using a variable separable technique, since boundaiy conditions 
have odd ordered derivatives in it, while the diflerential equation has 
even order derivatives Hence eneigy approach is used here, to obtain 
an approximate solution. The simplified energy criterion for this problem 
with lateral load p(*e,i/) is. 

Total potential, X[p = U + W — stationary value (4 14) 

Where U is the strain energy and W is the work done by external 
loads The total potential, lip, for a plate undergoing only bending can 
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be written as 


n. 


= u + w = 



^ d'^w 



+ 4£)(,() 


f 

[oT^j 



dx dtj 


= btationai y value 


(4 15) 


All the four edges of the lectangular plate consideied here are damped 
The boundaiy condihons along these edges are defined by the equations 


aUmq i = 0 and i = a jen all y 

w = 0 



do 


(416) 


aloiuj 1/ = 0 and y — b fo) all 'i 


iu — 


Oiu 

¥ 


0 

0 


(4 17) 


Rayleigh-Ritz method is used to solve Equation 4 15 A seiies approxi- 
mation foi the deflection of the plate «;()•,/;) can be defined as 




(418) 


The functions X,„(a;) and y;,(t/), that satisfy the boundary conditions 
4 16 and 4 17 are defined as follows 
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X,n(x) = ( ~ 

yitil/) = ( 2 /^ - ijy) V"~^^ (419) 


Substituting Equation 4 18 in 4 15 and applying the conditions given by 
Equations 4 16 and 4 17, we get M x N algebraic equations as 




+D 


12 


I dc2 ^ '"i) di/ d-c-* dy2 


, r. f<^dx,dx,„, fi>clY,ciy,,^\ 


dy dy 


ri» 

== / / V{'^,y)XmYnd'cdlJ 

Jq Jo 

ni = l,2, ,M and n = l,2, . 


(4.20) 


The system o[ M x N equations that we get after substituting Equation 
4 19 in 4 20 is of a form similar to 4 11 This can be transformed and 
solved for the unknowns afj as aheady described 


4 4 MIDPLANE SYMMETRIC LAMINATED PLATE 

We now consider plates with only bending-lwisling coupling The 
stiffness matrix teims present in the governing equation for w(r,!/) are 
— Dll, -^ 22 , ^ 12 / ^ 00 , Die and D 2 G Here again, we do not consider full 
equations along u° and since elements of [B\ matrix are identically 


zero. 
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4 4 1 Simply Supported Plate 


Foi midplane symmetric laminates, subjected to tiansverse static loads 
the Equation 4 3 leduce to the form 




. d^w 


+ 4-020 


d‘*w 

dxdy^ 


+ D 


S'* 


w 


22 


dy* 


= p(x,y) (4 21) 


The boundary conditions for the present case where all the toui edges 
aie simply supported aie 


along ® = 0 and x — a for all y 


tn = 0 


M, = 


^ d^w d^w d'^w 

~-^u-E-o- - 2 Dg(,o-;,_. ~ Di2 


dx^ 


dxdy 


= 0 


di/ 

(4 22) 


along y = 0 and y = b for all x 


w = 0 




d'^w 



^^'“axdy 

— Di 2 a 2 
Ol/ 

= 0 


(4.23) 


Closed form solution for Equation 4.21 is not feasible because of the 
presence of odd ordered derivatives in governing equahon and the 
boundary conditions. We formulate the problem in terms of the total 
potential ITp, and solve it by employmg Rayleigh-Ritz technique. Thus 
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the total potential Hj,, for the plate is. 


rr . TXT- 1 /■'’ /■" ^ /9‘iti V d^w d'w _ f d'w^ 

n„ - u + w - 2! i, V ^ "iv. 




£?'c2 

= stationary value 


(4 24) 


A senes solution of the form of Equation 4 18, that satisfies the boundary 
conditions specified by Equation 4 22 and 4 23 can be assumed here for 
solving Equation 4 24 Thus we get 


-\-D„ A„, y, dy+J^ X, <fa,y dy 

^D.J°X,X„,d. + ff^dy 

Jo Jo (hr cuir *^o aa; arc yo aiy ay 


H-2jDic 

+2D26 I [ 

\Jo 




0 (faS (i-c ‘^'/o dv^'-^'^k d% dx^ ‘^'io ^ 

• ^ <iA.^ /*rfy,<fy„ 

Xm-i—dv ' 


dF, ^ ^ ya ^ rfX,„ /'<' d% dYr, 

dv'''''Jo dy dy^ dx ^'^/o d^j^ dy 


dj/2 dy 


n a 

p['t,y)XmYndxdii 

. 

711 = 1,2, ,i\d and n = l,2, 


(4 25) 


M y N linear simultaneous equations obtained from Equation 4.25 is 
of the form specified by Equation 411. These set of equations can be 
solved for finding the statistics of the unknowns amn, if tbe statistics of 
the basic random variables 6f are known. The approximating functions 
given by Equation 4.8 are again used for solving the problem. 
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4 4 2 Clamped Plate 

Heie, again the goveimng equations and boundary conditions are given 
by Equations 4 21, 4 16 and 4 17 respectively We once again use 
the energy approach along with Rayleigh-Ritz technique satisfying the 
boundary condihons given by Equations 4 16 and 4 17 


4 5 NUMERICAL RESULTS 

Expciimental data regarding the distiibulion of basic composite material 
paiameleis is very limited in open literature. Analysis ot some of the 
available experimental data [2] show that the distribution of Poisson's 
ratio, Elastic modulus etc may be closely approximated by a Gaussian 
distiibulion For the present analysis, distiibution of all the landom 
variables are assumed to be Gaussian The knowledge of the distribution 
of the basic random variables is lequired for the Monte Carlo simulation 
approach. However, the peilurbation approach needs only the mean 
and SD of the vaiiablcs Table 41 gives the assumed mean values of 
all the primary random variables used. 


4 5 1 Validation of the Technique 


Validation of the proposed perturbation technique, for the static problem 
IS carried out for the case of a rectangular Graphite/Epoxy plate, with 
aspect ratio 1 5 by comparing the SD of deflection obtained with the 

centum, library 


I I r < AN PUR 

jmw 


Am ¥». A; i i 
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Table 4 1 Mean of Primary Random Variables [39] 


Random Variable 

Graphite/Epoxy 

Glass/Epoxy 

E\\ 

181 0 GPa 

53 78 GPa 

Eu 

10 3 GPa 

17 93 GPa 


0 28 

0 25 

Gi2 

7 17 GPa 

8 963 GPa 

<^23 

0 717 GPa 

0 89 GPa 


results of Monte Carlo simulation All the four edges aie assumed to be 
simply supported with each lamina having fibers oriented along X-axis 
The plate is loaded by a uiufoimly distiibuted load 

Results are computed for different sets of SDs for input RVs, to 
find the SD of deflection at the center of the plate The ratio aie 
assumed to be same foi all the foui primary RVs. 

Monte Carlo simulation solution is also obtained for the same prob- 
lem Two different random numbei generation techniques aie combined 
to generate the random samples with specified mean and SD for the 
system paramelcis lequned for the solution. This is done to reduce any 
bias m the random samples generated The total number of samples 
used m simulation is 15000 Figure 4 i shows the comparison of the 
proposed approach with the results obtamed from Monte Carlo simula- 
tion For the range of SD of input RVs considered the results from the 
present approximation comes very close to the Monte Carlo simulation 


results 




53 


The results show that the fust older perturbation technique adopted 
foi the piesent analysis is sufficient to give acemate lesiiUs, under 
the assumptions made foi the analysis of the pioblem foi the level of 
variations consideied in the basic random vniiables 

4 5 2 Simply Supported Rectangular Plate 

Results have been computed foi rectangular plates simply supported 
along all the foui edges The entire calculations aie earned out for 
plates ol constant thickness, h 

Fust some lesiills foi Graph ile/Epoxy laminates aie discussed To- 
wards the end some results for Glnss/Cpoxy laminates are also exam- 
ined Giaphile/Epoxy composites have latio aioiind 18 While foi 
Glass/Gpoxy composites it is of the order of 3 

Plots showing vaiialion of Mean of deflection and SD of deflection 
for one quailer of a plate are given m Figuie 4 2 The results are foi a 
[9()“] lamina, with aspect lalio of 1 0 All the input RVs were assumed 
to have SD equal lo 11) peicent of their mean values The nature of 
variation of SD acioss the plate is similar to the variation ol mean 
across the plate Points having higher magnitude of deflection have 
proportionately higher values of SD 

Figure 4 3 to 4 6 gives similar plots for the other basic RVs — the 
transverse modulus E 22 , the Poisson's ratio and the shear modulus 
(7 i 2 and longiludinal modulus En In all these cases SD of one of the 
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RVs IS changed and SD ol all othei RVs are kept zero 

All the foui plots show a lineai relationship between SD of the input 
RV and SD of deflection w at a point on the plate It is obseived that 
as we go towards the center of the plate the effect of input RV on the 
lesponse reduces Relation between distiibution of mean of deflection 
across the plate and SD of deflection across the plate is not linear The 
nature of this varies with the basic random variable that we considei 
The maximum eflect on the SD of deflection is due to change in Fn 
The effect of G 12 and i ?22 respectively, comes next Variation in SD of 
1^12 has the least effect on deflection Thus, if we consider a combination 
of effects of all the basic random variables, the net effect of them on 
the SD of deflection will be a complex combination of effects due to 
each one of them. The above observation implies that experiments for 
evaluation of Fu should be performed very precisely 

Figures 4.7 to 4.10 show the vaiiation of SD of coefficients a„i„ — m 
Equation 4.7 — with the change in SD of input RVs The results are 
for a [90“] lamina, with AR as a parameter The plots show a linear 
trend Increase in aspect ratio beyond a certain level doesn't seem to 
affect the results, As indicated by Figure 4 7, the first coefficient of the 
series an is affected most by the change in SD of input RVs 

Figures 4.11 gives the effect of SD of primary RVs on the SD of 
deflection ui, at four different locations for an aspect ratio of 0 5 Figures 
4.12 to 4.16 are for aspect ratios 1 0, 2.0, 3.0, 4 0 and 5.0 respectively 
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at different locations They also indicate a similai behavior The effect 
of increasing aspect ratio is to increase the SD of deflection This 
effect dominant at low AR, slows down with increase in AR Beyond a 
limit, the aspect ratio doesn't seem to have much effect on the SD of 
deflection 

Figure 4 17 shows the variation of SD of deflection w at the centei of 
the plate, with fibei orientation The results aie presented for different 
levels of SD of input RVs, ranging fiom 1 percent to 15 percent of the 
mean. At lower levels of input SD of RV, change in fibei orientation 
has practically no effect on the SD of deflection At highei levels of SD 
of input RVs the fiber orientation seems to affect the SD of deflection 
Fibei orientations which results in stiffei plates seem to have lower 
levels of SD of deflection. 

Variation of SD of deflection at the center of a plate, with the plate 
aspect ratio (a/b) is shown in Figure 4.18 foi [90“] stacking sequence 
Results are presented for different levels of SD of input RVs ranging 
from 1 percent of mean to 15 percent of mean The plate is assumed 
to be simply supported along all the four edges. As with the case of 
change in stacking sequence, here also change in aspect ratio doesn't 
have much effect on the SD of deflection at low levels of SD of input 
RVs. As the SD of input RVs is increased we can observe that there is 
significant influence of change in aspect ratio on the SD of deflection, 
at the initial stages of the curve But beyond a limit the result stabilizes 
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and further change in aspect ratio doesn't seem have any effect on the 
SD 

Change in SD of deflection with SD of input RVs is shown in Figiiies 
4 19 , 4 20 and 4 21 foi lammates having slacking sequences [0°] , [OO^/O^], 
and [0°/90°]^ respectively The results are foi simply supported square 
plates. These plots are similar to Figure 4 12 — the results for a [90°] 
square plate It may be noted here that all these plates have the same 
thickness and size. 

Figuie 4 22 and 4 23 show the plot of SD of deflection w against 
the SD of input RVs foi [45“] and {4r)‘'/0°], plates The plates compaied 
heie are square plates with identical thickness These plots aie almost 
identical — indicating that the middle 0° layers in the second plate 
does not have any noticeable influence on the result Results for 
[0°/45°/90“/— 45°]^ , [-45''/45°]g and [60°] laminates are given in Figures 
4 24 , 4 25 and 4 26 respectively The plots indicate a similarity in the 
behavior pattern in all the cases 

Now some typical results for a Glass/Epoxy composite are consid- 
ered. Table 4.1 lists the mean values of primary random variables used 
The variation of SD of deflection at different locations of the plate, 
with the SD of input RVs is given by Figures 4 27, 4.28 and 4 29, for 
aspect ratios 0.5, 1 0 and 2.0 respectively. The results are for rectan- 
gular laminates having [90°] stacking sequence, with all the four sides 
simply supported. The SD values are much larger than that for the 
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Graphite/Epoxy laminates having similar configuration, given by Figuies 
4 11, 4.12 and 4 13 This behavioi can be expected since the stiffness 
of Glass/Epoxy laminates are much less than that for Graphite/Epoxy 
laminates 

Figure 4 30 shows the sensitivity of SD of deflection at the midpoint 
of a [90°] lamina, to change in aspect ratio The lamina is assumed to 
be simply suppoited on all the four edges. Here the behavioi is similai 
to the the corresponding plot foi Graphite/Epoxy given by Figure 4 18 
except at the initial poilions of the cuives, coi responding to low aspect 
ratios At AR 1 0 the SD of deflection dips to a minimum value, aftei 
whicli jl gradually increases, with a tendency to readi a plateau at some 
high values of aspect ratio. Heie also, at low levels of SD of input RVs, 
change in AR do not have much influence on the SD of deflection 

4 5 3 Clamped Rectangjiar Plate 

Here we considei rectangular plates with all the foui sides clamped 
All the plates considered here are assumed to have the same thickness 
The results presented are for a Giaphite/Epoxy composite, specifications 
of which are given in Table 4 1 

The variation of mean of deflection and SD of deflection acioss the 
full plate, when all the edges are clamped is given by Figure 4.31. The 
results are for a [90°] square lamina, with SD of all the input RVs set to 
10 percent of their mean values. The nature of variation of the SD across 
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the plate is similar to that of mean Figures 4 32 , 4 33, 4 34 and 4 35 
show SD of deflection m against the change in SD of input RVs at various 
points across the plate, foi (O^’j , [90°j , [9070°]^ and [OV'^O'']^ laminates 
All the laminates under consideiation have the same thickness and 
platfoim dimensions The plots show identical tiends of vaiiation of 
deflection SD with the changes in the input SD The SD of deflection 
w, like the mean value is consideiably lower than the corresponding 
values for similai plates with simply suppoited boundary conditions as 
desciibed in the last section 

Results for [45°] and [45°/0°]j squaie plates are given in Figures 4 36 
and 4 37 respectively Unlike the results for similai plates with simply 
suppoited boundary conditions as discussed m the last section, heie 
the plots show difteient trends Here [45“/0°]g laminate shows a highei 
level of SD compared with the [45°] lamina Hence we can see that the 
middle 0° layers do have an influence on the SD of deflection for the 
present case 

Figures 4.38, 4 39 and 4.40 show the variation of SD of deflec- 
tion with the SD of input RVs, of laminates [60°] , [“46'^/45°], and 
[0°/45°/90°/— 45°]3 All of them show linear trends, with [~45°/45”)^ 
stacking sequence showing the minimum SD. This again shows the 
dominating influence of the outer plies, for lammates in the bending 
mode. 

The influence of change in fiber orientation on the SD of deflection 



59 


at the center ot a square plate is shown in Figuie 4 41 The results 
are piesented fot vaiious levels of SD of input RVs The trend here 
IS similai to that given in Figure 4 17 At low levels of SD of input 
RVs, the SD of deflection is not influenced much buy the change in 
fiber oiientation Fibei orientations which lesults in stiffer plates seem 
to have lowei levels of SD of deflection 


4 6 SUMMARY 

Application of the peituibation technique desciibed in Chapter 3 to the 
problem of deflection of leclangular plates with landomncss in basic 
mateiial properties is clemonslialed In this chaplei the formulation of 
the problem tor vaiioiis configuiations of boundary conditions, staclcing 
sequences, aspect ratios etc have been presented Results for various 
combinations of paiameteis involved have been presented The fust 
order perturbation foimulation and some of the associated assumptions 
made have been validated by an independent Monte Carlo simulation 
This analysis shows that the first oidei perturbation technique 
adopted for the displacement lesponse is sufficient to give accurate 
results foi the problem discussed In practice the overall effect is de- 
pendent on the extent of randomness exhibited by each one of the 
input random variables The maximum effect on SD of deflection is 
mainly due to change in E[i It is theiefoie necessary that experiments 
for evaluation of En should be performed veiy precisely to get ac- 
curate displacement icsponse estimation Most of the results show a 
linear trend for the lange of variations considered The Monte-Cailo 
simulation also exhibits a similar behavior Kence we may conclude 
that, apart fiom the effect of linear perturbation analysis, the response 
relations themselves arc linear in nature for the problem considered 
here. Fiber orientations which results in stiffer plates gives lower SD 
of deflection In general, at low levels of SD of input RVs, change in 
AR, fiber orientation etc doesn't have much influence on the SD of 
deflection The outer plies, in general have a dominating influence on 
SD of deflection 
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Figuie 41‘ Comparison of results from a Monte Carlo simulation with the 
present approximation Change of SD of deflection at the centre of the plate 
{0 5a, 0 5b), with SD of input RVs [0°] laminate, AR=1 5, all sides simply 
supported Graphite/Epoxy 



Figure 4 2: SD of input RVs 10 % of mean, [90“] laminate, AR=1 0, all sides 
simply supported. Distribution of mean of deflection and SD of deflection across 
of the plate. Graphite/Epoxy 
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Figure 4 3 Sensitivity of SD of deflection to SD of input RV I^ 2 > SD of all 
other input RVs kept zero [90“] laminate, AR=1 0, all sides simply supported 
Graphite/Epoxy 
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Figure 4.4: Sensitivity of SD of deflection to SD of input RV SD of al! other 
input RVs kept zero [90“], AR=1.0, all sides simply supported Graphite/Epoxy 



62 



cr» 

trt 



O QO O OS O 10 CIS O 20 

Input RV, SD y Mean 


Figuie 4 5 Sensitivity of SD of deflection to SD of input RV Gi^ SD of all 
other input RVs kept zero [90”] laminate, AR=1 0, all sides simply supported 
Graphite/ Epoxy 
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Figure 4.6: Sensitivity of SD of deflection to SD of input RV En SD of all 
other input RVs kept zero. [90°] laminate, AR=1 0, all sides simply supported 
Graphite/Epoxy 
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Figuie 4 7 Change in SD of syslem parameter an, with SD of input RVs [90°] 
laminate, all sides simply supported Graphite/Epoxy 
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Figure 4.8: Change in SD of system parameter aia, with SD of input RVs [90°] 
laminate, all sides simply supported Graphite/Epoxy 
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Figtiie 4 9 Change m SD of system parameter a^i, with SD of input RVs [90°] 
laminate, all sides simply supported Graphite/Epoxy 
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Figure 4. 10* Change in SD of system parameter aas, with SD of input RVs 
[90°] laminate, all sides simply supported Graphite/Epoxy 
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Lue 4 13 [90“j laminate, AR=2, all sides simply supported, Change in SD 

leflection with SD of input RVs Graph ite/Epoxy 
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Figure 4.14. [90“] laminate, AR=3, all sides simply supported. Change in SD 
of deflection with SD of input RVs Graph ite/Epoxy 
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.lie 4 15 [9(f] laminate, AR=4, all sides simply supported Change in SD 

leflection with SD of input RVs Graphite/Epoxy 
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Figure 4.16: [90°] laminate, AR=5, all sides simply supported. Change in SD 
of deflection with SD of input RVs. Graphite/Epoxy 
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Figure 4 19 [O'"] laminate, AR=:1, all sides simply supported Change in SD of 
deflection with SD of input RVs Graphite/Epoxy 
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Figure 4 20: [90“/0°]j laminate, AR=1, all sides simply supported Change in 
SD of deflection with SD of input RVs Graphite/Epoxy 





SD of W/ h SD of W/ h 


70 



O oo O OS O 1 O 0 15 

Input RV, SD/ Mean 


Lire 4.21 [0“/90“]g laminate, AR=1, all sides simply supported 

of deflection with SD of input RVs Graphite/Epoxy 
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Figuie 4 22: [45“] laminate, AR=1, all sides simply supported Change in SD 
of deflection with SD of input RVs Graphite/Epoxy 
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Liic 4 23 [4r)'’/0°], laminate, AR=1, all sides simply supported 

of deflection with SD of input RVs Graphite/Epoxy 
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Figure 4.24: [0745°/90'’/-45°]g laminate, AR=1, all sides simply supported 
Change in SD of deflection with SD of input RVs Graphite/Epoxy 
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Figure 4 25. [~4r)‘’/4r)°], laminate, AR=1, all sides simply supported Change 
in SD of deflection with SD of input RVs Graphite/ Epoxy 
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Figure 4.26’ [60“] laminate, AR=1, all sides simply supported Change in SD 
of deflection with SD of input RVs. Graphite/Epoxy 
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Figutu 4 27’ Test case Glass/Epoxy composite Change In SD of deflection, 
with SD of input RVs [90”) laminate, AR=0 5, all sides simply supported 
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Figure 4.28: Test case Glass/Epoxy composite Change in SD of deflection, 
with SD of input RVs (90°) laminate, AR=1.0, all sides simply supported 




Figuic 4 31 SD of input RVs 10 % of mean, [90°j laminate, AR=1 0, all sides 
clamped Distribution of mean of deflection and SD of deflection across the full 
plate Graphile/Epoxy 
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Figure 4.32. [0°] laminate, AR=1, all sides clamped Change in SD of deflection 
with SD of input RVs Graphite/ Epoxy 
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LU’c 4.33 [90“] laminate, AR=1, all sides clamped Change in SD of deflec- 
I with SD of input RVs Graphite/Epoxy 
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Figure 4 34. [90°/0“]„ laminate, AR=1, all sides clamped Change in SD of 
deflection with SD of input RVs Graphite/Epoxy 
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Figure 4 35. [()"/9()"],, laminate, AR=1, all sides clamped Change in SD of 
deflection with SD of input RVs Graphite/ Epoxy 
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Figure 4.36. [45“] laminate, AR=1, all sides clamped Change m SD of deflec- 
tion with SD of input RVs. Graphite/Epoxy. 
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Lire 4 37. [4r)'V0‘']^ laminate, AR=1, all sides clamped Change in SD of 
ection with SD of input RVs Graphite/Epoxy 
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Figure 4.38. [60“j iaminate, AR=i, all sides clamped Change in SD of deflec- 
tion with SD of input RVs Graphite/ Epoxy 
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Figure 4.41’ Vjiriation of SD of deflection lu at (0 5a. 0 5b). with fiber orien- 
tation for different SDs of input RVs. AR=:1, all edges of the plate clamped 
Graphito/Epoxy 



Chapter 5 


COMPOSITE PLATES — FREE VIBRATION 
ANALYSIS 

51 INTRODUCTION 

The* problem of reelangular plates undergoing fiee vibrations is con- 
sideied in this chapter. The classical laminated plate theory is used to 
obtain the governing equations. The randomness in basic material pa- 
lameteis is intioduced latei, using the techniques described in Chapter 
3 The gcometiy of the plate is defined by Figure 3.1. 

5 2 FORMULATION 

Consider a lectangular laminated composite plate undergoing free os- 
cillations The goveining equations of such a system in terms of dis- 
placements u", and w) and the stiffness matrix terms (Aj, A; and 
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Di'}) IS given by 140] 
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(5 3) 


Equtilions 5 1, 5 2 and 5 3 can be simplified based on stacking sequence. 


5.3 SPECIALLY ORTHOTROPIC LAMINATED PLATE 

In specially oilholropic laminates all the elements of [5] matux are 
identically zero. Other coupling terms like ^lioiAao and D\q,D 2 g are also 
absent. The only stiffness matrix terms present are Dui A 2 / D 12 and 
Dqq Again here we are concerned only with the transverse oscillations 
of the plate. Owing to these simplifications several terms drop out from 
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the goncial equalionb given above and icsuU in a simplified governing 
equation as 




71 , o/n , .-.71 s T. 0'K> 

-I 2(/5.. + + /’•^ = 0 (5 4) 




Closed foim solution foi equation 5 4 can be found only foi the case 
ol some specialized boundary conditions 


5 3 1 Simply Supported Plate 

lleie, we considei lectangulai composite plates with all the foui sides 
simply siippoited Even though closed foim solution can be found for 
the Lqualion 5 4 with this boundaiy condition, we tiy to attempt the 
solution using eneigy foimulalion. The simplified energy criterion for 
this problem assuming no laleial oi in-plane loads is 


Total pofcnlial, 11,, = U -\-T ~ stahonaiy value (5 5) 


Wheie II IS the sliain eneigy and T is the kmelic energy This for the 
piesenl case is: 
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(5.6) 
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Tlic boundniy condUions along the edges of the plate for the case 
when all the loin edges ol the plate aie simply suppoited aie defined 
by 


(Uoiuj I = 0 ai}d i = a for all y 

w 

M, 


0 

r, d^vj O^w 

I — “ — — ]D\*) 

0 (5 7) 


aUmy y 0 and y -* b for all 


w = 0 


My 


^ d'^w d^w 

0 (5 8) 


The Rayleigh-Rilz technique can be used to solve Equation 5 6 for the 
boundary conditions specihed by Equations 5 7 and 5 8 A series ap- 
proximation (oi the Iransveisc deflcclion iw('c,i/), satisfying the boundary 
conthlions can be defined as* 

Here X,n(‘i^) and Yniy) are eigen-funclions that satisfy the set of boundary 
conditions considered. Substituting Equation 5.9 in 5 6 and applying the 
conditions given by Equations 5.7 and 5.8 we get M xN homogeneous 
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algebraic equations 



The ap]UOximalin}5 funelion given by Equation 5 9 for deflection ^o{xyl/) 
salislying b.qualions 5 7 and 5 8 is assumed as* 




niTT'i, mry 

sin sm — — 

a b 


(511) 


Equation 5.11 can now be siibsliliiled in 510 to get the set oi M x N 
equations foi the present case. Equations 510 aie of the form of 
Equations 3 2 desciibed in Chapter 3 The proceduie described m 
Chapter 3 can now be applied here to mcoipoiate the randomness in 
basic material paramelcis and conveil the above equation to the form of 
Equation 3.1. These equations can be solved for the random eigenvalues 
using the random response w^{x,y), given by. 
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5 3 2 Clnmpecl Plate 


Ilcic all Iho loin edges o( the icclangulai plate consideied aie assumed 
to be clamped Ihe boundaiy condilioiis are defined by the ec|uations 


alou (] T — 0 and % = a for all y 


w ^ 0 


(513) 


(ilouy If - 0 and y = b for all x 


«) — 0 

P = 0 

dy 


(514) 


I’oi these boundaiy conditions closed foim solution of Equation 5 4 is 
dilficLilt. Again solution of the problem by energy methods is attempted 
The appioximatlng function foi — the tiansverse deflection of 

the plate — is again assumed as* 

«(■»,!/) = ( 515 ) 

The tu notions and y„(y), that satisfy the boundary conditions 

513 and 514 are defined again as follows: 


rM = 


(516) 


The governing equation for this case is also given by Equations 510 
Using the same procedure as carried out in the last section/ it can 
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be shown that, the lesulhng equation is of the foim ot Equation 3 2 
fioin Cliaplei 3 The landom appioximating function toi deflection is 
assumed to have the foim 

(517) 

Using this. Equation 5 10 can be solved as desciibed in the last section, 
foi the random eigenvalues, resulting fiom the randomness m basic 
maleiial paiamelers. 


5 4 MIDPLANE SYMMETRIC LAMINATED PLATE 

Bonding-twisting aie the only coupling terms that are present in mid- 
planc symmetric angle-ply laminates. !?(,, T>) 2 / Z>(,g/ T'lo and 

aie the stiffness matrix lei ms piesenl in the governing equation In 
Equation 5 3 we arc concerned only with the w terms, further [B] inatiix 
terms aic zeio. This results in a much simplified governing equation 


5 41 Simply Supported Plate 

For midplane symmeliic laminates Equation 5.3 takes the following 
torm: 


Dn 


d^w 

dfc'' 


ID 

j~, d'^w d^V) 


= 0 (5 16) 
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’IhL- boundaiy loiidilions aii> ajviin given by Equalions 4 22 and 4 23 
The appioximaliug funelion (oi w{i, ij) can again be assumed to be of 
llie foim of I'.ejualion 5 U Again, closed foun solution of this Equation 
5 18 IS not (easiblc If, solution of the problem is attempted using eneigy 
foimulation. Die total potential encigy II,, foi the plate is 
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(5 19) 


The appioximating lunction given by Equation 5.11 along with Equation 
5 19 IS used to lind the governing equations — using Rayleigh-Ritz 
lechnK|ue Heme we gel 
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This ‘.cl ol A/ X N homo^’cncous sitnuUnncoiis equations can be expiessecl 
jn Ihc foun of J'qiiiUioii 3 8 given in Chti|itci 3 The cleiiveci equations 
foi finding the mean and vaiianre of the eigen values can be expiessed 
in the ioim of ['.quations 3 10 and 3 11 Foi a non-luvial solution of 
these equations, eoellicienls n.j cannot be zeio Hence lest of the teims 
in these system of equations must go to 7cio This condition will give a 
set ol equations that loiin an eigenvalue pioblein When solved, these 
sets ol equations yield men and vaiiance of noimal mode frequencies 

5 4 2 Clain[)LHl Plate 

The tlamped boundaiy eondilions aie specified by the Equations 5 13 
and 5.14 The appioximaling function assumed is again given by Equa- 
tions 5 15 and 5 U> I he governing equations foi the conventional 
pioblem IS ome again }',iven by Equations 5 20 The landomness is 
incoipoialed in these equations as desciibed in Chaplei 3 


5 5 NUMERICAL RESULTS 

Results have been obtained foi some typical pioblems to study the effect 
of vaiiation m nialeual propeily chaiactciistics on the natuial fiequencies 
of K’llangulai eoinposile plale.s All the plate configuiations discussed 
aie assumed to have the same total thickness Foi the pioblems studied, 
the basic inalcuial ptopedies like the longitudinal modulus, transvcise 
modulus, sheai modulus etc. are considered to be random vanables. 
The landomness in maleiial piopeilies is handled using the peiturbation 
Icchnicpio discussed in Chapter 3 and earlier sections in this chapter 
The assumed mean values of the landom vanables used aie given in 
Table 4 L of Chaptei 4 
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5 5 1 Vjlidalioti of the Tecliniqup 

Vdlid^Uioii ol (lie perliiibolion lec!inu|iie u^ed, for the piesenf ptoblem 
was done by compaiiiifj il with Ihe icbults fiom a Monle Catlo simulalion 
study The basic random vniuibles aie assumed lo be Gaussian for Ihis 
puipose Since solving hqualion 5.11 using Monte Cailo simulation is 
compulnlionnlly expensive, an one leim appioximation of the equation 
was used for validation ( Results obtained by solving the full set of 
equation is given la lei ) 

l’i}',uie 5 I shows the results liom the Monte Carlo simulation, com- 
pared with the present approximation for a [90“j square Graphitc/Epoxy 
plate with all the lour sides simply suppoited The ratio are 
assumed to be the same lor all the foui piimary random variables. 
Random samples ol the piimary RVs with specified chaiacteristics weie 
generated lor simulation, as described in Chapter 4. Here the SD of 
fundamental and higher natural frequencies are plotted against the 
change in SO of input KV.s. As indicated by the plots, for the range of 
SD of input RVs considcicd both the results are quite close. One may, 
thereloie, conclude that the present approximation employing first or- 
der peiturbalion gives sufficiently acciuale results for the free vibration 
pioblem under consideration 
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5 5 2 Sunply Suppodocl Rcclcjngular Plate 

In this sc'Llion we I’ousulei fiec vibiation of thin reclangulai composite 
plnloh himply suppiuletl along all the loiii' edges Except when indicated, 
all the input KVs aio assumed to have the same SD All the plates undei 
consideiation aio assumed to have the same thickness h Most of the 
lesults aie ten Ciraphtte/Hpoxy plates Towards the end some lesults foi 
plates made ot C^lass/Lipoxy composite aie also presented This is done 
to jnvestip,ali’ the ellecl of modulai lalio since the ratio of is more 
foi Cliaphile/lipoxy malenai, in compaiison with Glass/Epoxy 

The mean values ol natural hequencies toi a [90“] lamina simply 
buppoited along all the edges is given in Table 51, foi various aspect 
ratios lable 5 2 gives mean of natural Irequencies for a square laminate, 
with vniious stacking sequences Figure 5 2 gives the plot of variation 
of SD of natiual fiequencies of a [90“) square lamina with variation 

o 

in SD of longitudinal modulus l!?n SD of all other RVs are assumed to 
be zero. F'iguies 5.3 , ^5.4 and 5.5 show similar results for the other RVs 
— llie Iransvenw modulus, the Poisson's ralio and Gi 2 the sheai 
modulus The SD of natural fiequcncics arc affected most by the change 
in SD of the longitudinal modulus E\\ The fundamental frequency is 
most affected by changes in JSn, E'n and i>u • We can see that the 
influence of SD of input RVs on SD of natural frequencies diminishes, 
as we consider higher and higlier modes for the above three cases. We 
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can luilluM obseivo a flallcning Ueiui foi thib effect, as indicated by 
I'lgutes 5.6 , 5 7 and 5 8. In tho case of G\i, (he fmidaiuental liequency 
IS least allccled by (lie change in SI9 of inpiil RV The highei modes 
aie piogiessively a ducted moie by the change in SD of Gia Here also 
theie IS a flattening liend loi this cited as indicated by Figiue 5 9 It 
we consider the magnitude ot edecl on the SD of natuial fiequencies 
by the SI) each ol the input RVs — En, G 12 , E 12 and respectively 
show piogicHsively less and less inlluence on the SD of the natural 
fiequencies 

I’lgiiie 5 U) shows the influence of simultaneous change m the SD of 
input RVs on the SI) ol the naliiial fiequencies. Here all the basic RVs 
aie assumed to have the same level of SD The results are for a (90“] 
lamina with aspect ratio 0.5 All the four sides aie assumed to be simply 
supported. The fust natural ticquency is affected least. Higher and 
highei modes are piogiessivcly aflected more by the change in SD of 
input RVs. Tigiires 5.ii , 5.12 , 5.13 , 5.14 and 5.15 give corresponding 
results lor aspect ratios 1.0, 2.0, 3.0, 4.0 and 5.0 respectively. In all 
these ca.ses, except tor the aspect ratio and SD of input RVs all othei 
parameters lomain the same. 

The influence of change in aspect ratio on the SD of natural frequen- 
cies IS shown by Figure 5.16 . Here all the input RVs are assumed to 
have a SD equal to 10% of their mean value. The results are for a [90“] 
lamina. The nature of variation of the SD of frequencies at low aspect 
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latios it> stiongly d0pendent on the mode we are considering But the 
lesuU seems to leach plateau as we appioach higher aspect rahos At 
aspect ratios above 2 0, the SD of first mode frequency is the largest, 
with gradually decreasuig values for the hrgher modes compared to the 
first mode shown in the figure At low aspect ratios this pattern is not 
present 

Figure 5 17 illustrates the influence of change in aspect ratio on the 
SD of fundamental frequency at differeirt SDs ot input RVs, with all of 
them having the same SD When the SD of input RVs is low the change 
in aspect ratio doesn't seem to affect the SD of natural frequencies 
As the SD of input RVs is increased, the influence of change m aspect 
ratio incieascs. The change in aspect ratio doesn't seem to affect the 
SD of 01 “ beyond a limit, depending on the SD of input RVs 

The influence of fiber orientation on the SD of fundamental frequency 
IS illustrated in Figure 5.18 — foi a squaie plate The influence of 
variation in SD of input RVs is also shown Here also for low levels of 
SD of input RVs, change in fiber orientation does not have much effect 
orr the SD of natural frequencies. For higher levels of SD of input RVs 
the SD of u)'^ uses to a maximum value at 45° fiber orientation and then 
again drops down. 

Results for the case of a [9070°]. square plate is shown m Figure 5.19 
for several frequencies. The results shown here are almost identical to 
result as given by Figure 5 11 for a [90°] lamina with the same dimensions 
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and thickness Figures 5 20 and 5 21 gives similar results for [0790°], 
and [0°] laminates 

Figures 5 22 and 5 23 shows the results for [45°] and [45°/0°], lami- 
nates The plots foi different modes are closely packed together m both 
the cases The 0° layers in the second case doesnT seem to have much 
mfluence on the output SD 

Results for [60°] , [-45745°], and [0°/45°/90°/-45°], laminate con- 
figurations are shown in Figuies 5 24 , 5 25 and 5 26 respechvely for a 
square plate In Figures 5 25 and 5 26 only the fundamental frequency 
is given, as other hequencies are very close. 

Some results for a typical Glass/Epoxy composite are now presented 
Table 5.3 gives mean values of natural frequencies for a [90°] lamina, 
simply supported along all the four edges Tlie result is given for 
different aspect ratios Figures 5 27, 5.28 and 5 29 give the variation 
of SD of natural frequencies with the SD of mput RVs for aspect 
ratios 0.5, 1.0 and 2.0 respectively The results are for a [90°] lamina, 
simply supported on all the four edges. The magnitude of SD of the 
natural frequencies are comparable to those for a similar Graphite/Epoxy 
configuration given by Figures 5 10, 5.11 and 5.12, except for the order 
of rate of change for the SD of the frequencies. 

Influence of the change m aspect ratio of the plate on the SD of 
natural frequencies is shown by Figure 5 30, for different levels of SD 
of input RVs. The result is for a [90°] lamina with all the sides simply 
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suppoiled Al low lev6ls of SD of input RVs the change in aspect ratio 
has piaclically no effect on the SD of first natural frequency At highei 
levels of SD of input RVs the response SD dips to the lowest value at 
AR 1 0, to to again inciease to higher values All the response curves 
ultimately reach a plateau and further increase in aspect ratio doesn't 
seem to have much effect on the SD of fiist natural frequency 

5 5 3 Clamped Rectangular Plate 

Heic we consider rectangular plates with all the four sides clamped 
All the plates aie assumed to have the same thickness 

Even though the numerical values are diffeient, the trend here also is 
similar to the simply suppoited boundary condition examples discussed 
in the last section The mean values of fundamental frequencies for 
square laminates with various stacking sequences are given in Table 
5 4 The results for variation in SD of with SD of input RV are 
shown in Figures 5.31, 5.32, 5.33, 5 34, 5 35 and 5 36 for the first natural 
fiequcncy for diffeient stacking sequences for a plate with aspect ratio 
1.0 In all these cases the higher modes are closely packed along with 
the fundamental frequency. All of them show linear behavior Figure 
5.37 illustrates the influence of variation m fiber orientation on the SD 
of the fundamental frequency for a square plate, for different SDs of 
input RVs. 

The rate of change as well as the order of change of the the natural 
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frequencies with change in SD of input RVs is dependent on the the 
stacking sequence and the boundary conditions of the plate i e for 
different configurations they vary differently This can be important 
when the plates aie subjected to an excitahon which approaches the 
range of variation of the natural frequencies 


56 SUMMARY 

Application of the perturbation technique described in Chapter 3, to 
problems free vibrahon of lectangular plates has been discussed For- 
mulation and results aie presented for problems with clamped and 
simply supported boundary conditions, having various combmahons of 
staclong sequence and aspect ratio The method has been validated for 
the piesent class of problems by an independent Monte Carlo simula- 
tion, which shows excellent agreement with the results from the present 
technique 

The perturbation approximation gives suteciently accurate results for 
the problem under consideiation In general the SD of natural frequen- 
cies are affected most by the change m En The nature of influence of 
the input RVs are strongly dependent on the mode of vibration and the 
extent of randomness of the variable under consideration The nature 
of variation of SD of normal mode fiequencies is strongly dependent on 
the mode under consideration at low aspect ratios Change m aspect 
ratio beyond a limit doesn't seem to affect the SD of natural fiequencies 
At low levels of SD of input RVs change m aspect ratio, fiber orientation 
etc doesn't seem to have much effect on the SD of natural frequencies 
SD of natural frequencies have a maximum value for a laminate with 
45° fiber orientation The outer plies of laminates have a dominating 
influence on the SD of natural frequencies In general the lesults show 
a linear behavior foi the change in SD of natural frequencies with the 
SD of the input RVs The rate of change and order of change of natural 
frequencies with change m SD of input RVs is strongly dependent on 
the stacking sequence and the boundary conditions. 
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Table 5 1 Mean of natural frequencies for different aspect ratios [90*^] 
lamina, all sides simply supported Graphite/ Epoxy 


AR 


2nd 

3 rd 



0 15044E-02 

0 32505E-02 

0 61606E-02 

0 10235E-01 


0 10679E-02 

0 15044E-02 

0 22319E-02 

0 32505E-02 


0 95878E-03 

0 10679E-02 

0 12498E-02 

0 15044E-02 

30 

0 93857E-03 

0 98707E-03 

0 10679E-02 

011811E-02 

40 

0 93149E-03 

0.95878E-03 

0 10042E-02 

0 10679E-02 

50 

0 92822E-03 

0 94568E-03 

0 97478E-03 

0 10155E-02 


Table 5 2. Mean of natural frequencies uP', for different stackling sequences 
Square plate, all sides simply supported Graphite/Epoxy 


Stacking Sequence 


2iu/ 

3rd 


190VO-1. 

0 959E-03 

0.139E-02 

0 212E-02 . 

0 3I4E-02 

(0790”|, 

0 306E-03 

0 743E-03 

0147E-02 ' 

0.248E-02 

[45“| 

0 123E-02 

0 309E+00 

0 393E+00 

0 411E+00 

[60»1 

0 130E-02 

0 361E+00 

0 402E+00 1 

0 431E+00 

14570"], 

0 llOE-02 

0 277E-)-00 

0 352E+00 : 

0 369E+00 

(-46«/45"], 

0 123E-02 

0.318E+00 

0.393E+00 

0 410E4-00 

(0°/46°/90"/ - 45°], 

0 616E-03 

0177E+00 

0 190E+00 

0 210E+00 
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Table 5 3 Mean of natural frequencies for different aspect ratios [90°] 
lamina, all sides simply supported Glass/Epoxy 



AR=0 5 

AR=1 0 

AR=2.0 


0 10062E-02 

0 46055E-03 

0 32412E-03 


0 31890E-02 

0 10062E-02 

0 46055E-03 

3.d 

0 68269E-02 

0 19157E-02 

0 66792E-03 

4** 

0 73687E-02 

0 31890E-02 

0 10062E-02 


Table 5 4 Mean of fundamental frequency itP", for different stackling sequences 
Square plate, all sides clamped Graphite/ Epoxy 


Stacking Sequence 

0.2 

[90“] 

0 118E-04 

190”/0”], 

0 118E-04 

(0V90"1. 

0.118E-04 

[45"] 

0 378E-04 

[60") 

0 334E-04 

[ 45701 , 

0 346E-04 

[- 45745 °], 

0 662E-05 

[0°/45°/90°/ - 45°}, 

0 193E-04 
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Figure 5 !• Normalized SD of natural frequencies plotted against SD of 
input RVs The present approximation compared with results from a Monte 
Carlo simulation [90°] laminate, AR=1 0 Graphite/Epoxy 
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Figure 5 2. [90°] laminate, AR=1, all sides simply supported, sensitivity of SD 
of normalized natural frequencies to SD of input RV E\[. SD of all other 
input RVs kept zero Graph ite/Epoxy 
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Figure 5 3 [90°] laminate, AR=1, all sides simply supported, sensitivity of SD 
of normalized natural frequencies to SD of input RV E 22 SD of all other 
input RVs kept zero Graphite/Epoxy 
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Figure 5 4: [90°] laminate, AR=1, all sides simply supported, sensitivity of SD 
of normalized natural frequencies to SD of input RV SD of all other 
input RVs kept zero Graphite/Epoxy 
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Figiiie 5.11 [90°] laminate, AR=1, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/Epoxy 
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Figure 5 12, [90''] laminate, AR=2, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/ Epoxy 
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Figure 5 13 [90°] laminate, AR=3, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/Epoxy 
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Figure 5 14 [90°] laminate, AR=4, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/ Epoxy 
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Figure 5 15 [90"^] lamina, AR=5, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/ Epoxy 
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Figure 5.16: [90°] laminate, AR=1, all sides simply supported Variation of SD 
of normalized natural frequencies with change in aspect ratio, for different 
modes SD of input RVs 10 % of mean Graphite/Epoxy 
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Figure 5 17. [90°] laminate, all sides simply supported Variation of SD of 
normalized fundamental frequency with change in aspect ratio, for different 
SDs of input RVs Graph ite/Epoxy 
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Figure 5 18 AR=1, all sides simply supported Variation of SD of normalized 
first natural frequency with change fiber orientation, for different SDs of input 
RVs Graph Ite/Epoxy 
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Figure 5 19 [QO^'/O"]^ laminate, AR=1, all sides sinnply supported Variation of 
SD of normalized natural frequencies with SD of input RVs Graphite/Epoxy 
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Figure 5.20, [0''/90'’]^ laminate, AR=1, all sides simply supported Variation of 
SD of normalized natural frequencies with SD of input RVs Graph ite/Epoxy 



109 



CO 


cvj •' 



O OO O 05 O 1 O O 1 S O 20 

Input RV, SD / Mean 


Figuie 5 21 [0°] laminate, AR=1, all sides simply supported Variation of SD 

of normalized natural frequencies oj^ with SD of input RVs Graphite/ Epoxy 
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Figure 5-22 [45°] laminate, AR=1, all sides sinnply supported Variation of SD 
of normalized natural frequencies uP' with SD of input RVs Graphite/Epoxy 
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Figure 5.23 [45°/0°|j laminate, AR=1, all sides simply supported Variation of 
SD of normalized natural frequencies uj'^ with SD of input RVs Graphite/Epoxy 



Figure 5,24’ [60°j laminate, AR=1, all sides simply supported Variation of SD 
of normalized natural frequencies with SD of input RVs Graphite/ Epoxy 
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Figure 5,25 [~45°/45®]j laminate, AR=1, all sides simply supported Vari- 

ation of SD of normalized fundamental frequency oj"^ with SD of input RVs 
Graphite/Epoxy 



Figure 5 26* [0°/45'’/90'’/-45'’j^ laminate, AR=1, all sides simply supported 
Variation of SD of normalized fundamental frequency with SD of input RVs 
Graphite/Epoxy 
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re 5.29 Test case glass/epoxy composite [90°] laminate, AR=2 0, all 
simply supported Variation of SD of normalized natural frequencies 
SD of input RVs 



OO 10 20 30 40 so 60 

aspect ratro, a/b 


Figure 5,30, Test case glass/epoxy composite. [90°] laminate, all sides simply 
supported Variation of SD of normalized fundamental frequency oj^ with change 
in aspect ratio, for different SDs of input RVs 
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Figure 5 31 [45°] laminate, AR=1, all sides clamped Variation of SD of 

normalized fundamental frequency with SD of input RVs Graphite/Epoxy 
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Figure 5 32 (60"] laminate, AR=1, all sides clamped Variation of SD of 

normalized fundamental frequency with SD of input RVs Graphite/Epoxy 
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Figure 5 33 [45“/0°]j laminate, AR=1, all sides clamped Variation of SD of 

normalized fundamental frequency with SD of input RVs Graphite/Epoxy 
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Figure 5.34. [-45“/45“]j laminate, AR=1, all sides clamped Variation of SD 
of normalized fundamental frequency u)^ with SD of input RVs Graphite/Epoxy 
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Figure 5 35 [90°] laminate, AR=:1, all sides clamped Variation of SD of 

normalized fundamental frequency ux^ with SD of input RVs Graphite/Epoxy 




Figure 5 36 [0°/45°/90“/'-45°]j, laminate, AR=1, all sides clamped Vari- 

ation of SD of normalized fundamental frequency with SD of input RVs 
Graphite/Epoxy 
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Figure 5 37 AR=1, all sides clamped Variation of SD of normalized first 

natural frequency with change fiber onentationi for different SDs of input 
RVs Graph ite/Epoxy 



Chapter 6 


COMPOSITE PLATES — BUCKLING ANALYSIS 

6 1 INTRODUCTION 

The aim of this chapter is to investigate the effect of material propeity 
randomness on the initial buckling characteiistics of rectangular com- 
posite plates The geometry of the plate and the co-ordinate system 
used are given in Figure 3 1 The classical laminated plate theory is 
used for the basic formulahon, to arrive at the governing equations for 
the problem 

6 2 FORMULATION 

Consider a rectangular lammated composite plate subjected to uniform 
inplane compressive and shear loads The plates considered are sup- 
ported along the all the four edges The governing equations for such a 
system in terms of displacements v° and w) and the stiffness matrix 
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terms (>1,^, and Dij) are given by [40] These are 


+ 2 A^q-— ~ - + -4.66 -?r-T- + Ai6 


dx^- 


dx"^ 


dxdy '’“92/2 

+ (-4i 2 + -Aee) oia.. + -^^e-PTr ~ 


dxdy 


^2/2 


dx^ 


. _ 9^1/; , _ , 92to d^w 

~ 3-5i 6 Q ,n. “ (-Bi2 + 2566) q..q. -, ~ ^26' 


dx^dy 


dxdy"^ 


dif 


= 0 


( 61 ) 


92 92ti® 92it° 92?;° 

^i6'^~r + (^12 + -^ee) Q_a.: + ^26"S75“ + ^ee 


dx 


dx'^ 


dxdy ' 9 j /2 

n A A -r, / Ti « -n \ 

+ 3^26 Q n +Vl22-^ 5 i6-^-7 — (5 12 + 2566) „ 2 a 

dxdy dy^ dx^ dx^dy 

d^w d^w 


d'^w _ 9'*to -,_ _ . d'^w _ d'^w 

Al-^^ + 45l6flj3a:.’ + 2(5 i2 + 2566) Q-.oo o + 4526 


dx^dy 


dx^dxj^ 


d^w 


93u° 


+^“■§7" 


“ 35i6 


92 gO 

dx^dy 


~~ (5i2 + 2566) 


dxdy^ 

92^° 


dxdy^ 


d^u^ 92^0 92^0 

_ S5 .^f--B — = iV — +27\r + r6 31 

9a;9y2 ^22 ^^3 ^^2 + 2 i\r^y + Ny (6 3 ) 


Equations 6 1, 6 2 and 6 3 can be simplified based on stacking sequence 


6 3 SPECIALLY ORTHOTROPIC LAMINATED PLATE 

For the case of specially orthotropic laminates the entire [5] matrix 
elements reduce to zero Other coupling terms like -Aic = A26 = 0 and 
Djg = D 26 = 0 also vanish So the only stiffness matrix terms present 
are Du, D 22 , D 12 and 566 Hence^ for the present case, the general 
equations given by Equations 61,62 and 6.3 are reduced to the form 
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Ai?? + 2{A2 + 2£)66)- 


dx^ 

d^w 


dx^-dy- 


+ -0: 


22 


d^w 


ly 


dxdy '*' ^ dy'^ 


(6 4) 


We consider the case of uniform compressive load applied along two 
opposite edges, in the X-direction Hence Equation 6 4 further reduces 
to the form 


i5u^+2(A.+2D„)^ 


+ D 


22 


d^w 

w 


= K 




w 


dx"^ 


(6 5) 


6 3 1 Simply Supported Plate 

The boundary conditions for the case when all the four sides of the plate 
are simply supported are defined by the following set of equations 


along cc = 0 and x = a for all y 


w ~ 0 




d^w 


~Du « ■ W ~ D\2 


d^w 

dy^ 


= 0 


(6 6 ) 


along y = 0 and y = b for all x 


w ~ 0 

d'^w d^w 

My = 


0 


(6.7) 
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A senes approximation for the deflection of the plate, w{x, y) is written 
as* 


w{x,y) = (6 8 ) 

where and Ynip) are functions that satishes the set of boundai/ 

conditions considered given by Equations 6 6 and 6 7 In terms of beam 
characteristic functions, Equation 6 8 is rewritten as 




rmrx niry 
mn Sin sin 


If Equation 6 9 is substituted in the governing Equation 6 5 and simpli- 
fied, the expression for the initial buckling load is obtained in the form 
[40] 

2 

No = [Dnm^ + 2{Di2 + 2DeB)m‘^n'^R'^ + D22n^R‘^] (6 10 ) 


Equation 6 10 reduces to the critical buckling load expression when 
n = 1 That is 


2 

Na = [x>iini‘' + 2{Di 2 + 2Dm)m^R'^ + D22R^] (6 11) 


Equation 6.11 is of the form of Equation 3 2 described in Chapter 3 At 
this stage we introduce the effect of randomness in material properties 
on the system response The procedure described m Chapter 3 is now 
applied here to incorporate the randomness in basic material parameters 
and convert the above equation to the form of Equation 3.1 These 
equations are then solved for the random critical buckling load N^, 
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using the random response y), given by 


w‘ 






miTx UTvy 
sm sin — — 


(612) 


The mean and variance of the random critical buckling load iVj can 
be found, provided the statistics of the basic random variables are 
known 


6 3 2 Clamped Plate 

For the case of a plate clamped on aU the edges it is difficult to solve 
the governing Equation 6 5 in closed form So the energy method 
IS used here to formulate the problem and to find an approximate 
solution Simplified energy criterion for the present problem, with 
inplane compressive loads is given by 

Total 'potential Tip — U -VW = stationary value (6 13) 


Where U is the strain energy, W is the work done by the external loads 
For the present case total potential 11^ is given by 

d'^in d^w 


1 /'*' /■" ^ f d'^w\ ^ 

n, - U+W = 4-213. 


+-D22 




+ 4i)( 


66 


/ d^w Y 

\dxdy ) 




Y 

^dx^j 


dx"^ 


dxdy 


= stationary value 


(614) 
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The boundary conditions for the case when all the sides are clamped 
are given by 


along a? = 0 and x = a for all y 

w = 0 

^ = 0 (615) 

OX 


along y = 0 and y ~ b for all x 


w = 0 
du) 


(616) 


We employ Rayleigh-Ritz technique to sole the Equation 6 14 for the 
boundary conditions specified by Equations 615 and 616 A senes 
approximating function for w{x,y) satisfying the boundary conditions 
given by Equations 6 15 and 616 can be defined as 




(6 17) 


with and Yn(y) that satisfy the boundary conditions given by 

= (a;^ — 

y ,.( y ) = (618) 


Substituting Equation 6 17 in Equation 6 14 and applying the conditions 
given by Equations 6 15 and 6.16 we get the following set of M x iV 


equations 
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rl’ 

d^X, “'■ 


i 




iiy2 


rdx,dx^^ r^dy^dv^^ 

■hD22 X.X^dx -j~f-j-^dy+4I>^Q y-dx , ■• <Xy 

JO Jo ay-^ dy-^ Jo dx dx Jo dy dy 

m — 1,2, ,M and n = l,2, ,N (619) 


Equahon 6 19 can once again be expressed in the form of Equation 3 2 
described in Chapter 3 The effect of randomness in primary random 
variables can be introduced m Equation 6.19 using the techniques de- 
scribed in Chaptei 3 The characteristics of random critical buckling 
load can now be found, if the statistics of input random variables 
bf are known 


6 4 MIDPLANE SYMMETRIC LAMINATED PLATE 

Midplane symmetric laminates have the bending-twisting coupling 
terms present in the governing equations The coupling stiffness matrix 
terns [S] are identically zero. The stiffness matrix terms Du, D 22 r 
Dee/ Die and D 26 are present in the governing equation 



125 


6 4 1 Simply Supported Plate 


For midplane symmetric laminated plates the Equation 6 3 reduces to 
the following form 


_ d'^w _ 


dx^dy 

^ ^ d'^w ^ d*w 
+ 4 £> 26 o ' o -i + -^22 


d^W 

■P 2(Z?i2 + 2Z166): 


dx‘^dy^ 


N. 


dxdy^ 

(9^10 

dx^ 




(6 20) 


The boundary conditions when all the four sides are simply supported 
are given by Equations 4 22 and 4 23 Since closed form solution of 
Equation 6 20 is difficult for this set of boundary condihons, we attempt 
a solution of the problem by energy method Here the total potential 


lip can be expressed in the following form 


= uj: 


f d^w\^ 

^lil I + 2Z) 




J 


12' 


'^D22 

-h4I?66 


d 2 ' 

-w 


dxf~) 

( d~w \ 


+ 4A6 
2 


+ iV, 


dxdy J 
= stationary value 


d'^w d^vj 
dx- dxdy 

m) 


+ 4£>; 


26 


dx^ di/ 

d'^w d^w 
dy"^ dxdy 


dxdy 


(6 21 ) 


The approximabng function for iu{x,y) satisfying the boundaiy condi- 
tions can be assumed again m the form of Equation 6 9 From these 
equations applying Rayleigh-Ritz method we get a set of M x iV linear 


equations as* 
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(PXid^Xrji ^ J 

E,=l E,=1 |-C^11 j[, ^^2 ^i^ndy 
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dx^ dx^ 

■ r ^ d ^ x ,^ rK .^ yn , r . 

/ -^m-T^rfa: / iS“7T^2/+ / X^~r~^dx / Y^-^dy 
JQ dx^ JQ dy^ Jo dx^ Jo dy^ 


^Dj‘ X.X.& / r /‘ 

./o JO dy^ dy^ Jo dx dx Jo dy dy 


+2£>i6 
+ 2£>28 


rdXd^x^_^^ 


da; da;2 


fx„^dx fJ^^dyd- fx,^dx J 

Jo dx Jo dy d'lr Jo dx Jo 


Jo 


dy dxp' 

m = l,2, ,M and n = l,2, ,i\r 


f‘> d% dY^ 
dy^ dy 


dy 


(6 22 ) 


Equation 6 22 is of the form of Equation 3 2 discussed m Chapter 3 
The effect due to randomness in basic material parameters can now be 
introduced by assuming an approximating function as given below for 
the random process for deflection w^{x,y) 


From Equation 6 22 and 6 23 employmg the technique described in 
Chapter 3 statistics of the random critical buckling load can be 
found, provided the statistics of the primary random variables are 
known 


6 4 2 Clamped Plate 

Here again the energy formulation along with the Rayleigh-Ritz tech- 
nique IS used to solve the problem The total potential energy equation 
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and the derived governing equation are given by Equations 6 21 and 
6.22 respectively The boundary conditions for the case where all the 
four sides of a rectangular plate are clamped are given by Equations 
6 15 and 6 16 Equations 6 17 and 6 18 give the approximating functions 
satisfying these boundary condihons A senes appioximating function 
like Equation 6 23 is used to incorporate the randomness owing to the 
basic material characteristics The characteristics of the random critical 
buckling load process is found using the procedure described in 
Chapter 3 

6 5 NUMERICAL RESULTS 

The effect of randomness in mateiial parameters on the buckling load 
and the mode shape is described in this section We limit our discussion 
to rectangular plates Thickness of the plates are taken to be the 
same irrespective of the change in laminate configurations. The mean 
values for the basic material parameters , E 22 » f?i 2 and uu for 
Graphite/Epoxy composite are as given in Table 4.1 The randomness 
in these basic material parameters is mcorporated by the perturbation 
technique described m Chapter 3 Bulk of the results presented are 
for a typical Graphite/Epoxy composite. Some results are given for a 
Glass/Epoxy composite system also 
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6 5 1 Validation of the Technique 

Validation of the present technique applied to the plate buckling problem 
IS done by comparing the results obtained using the present technique 
with the results from an independent Monte Carlo simulation The 
problem of a [90“] rectangular lamina with aU the four sides simply 
supported, subjected to a uniformly distributed compressive load from 
two opposite X-edges is considered for this The results are computed 
for a Graphite/Epoxy system In Figure 61 SD of normalized buckling 
load Ncr {Ncr = is plotted against the aspect ratio for various SD 
of input RVs ranging from 1 percent to 15 percent of mean Results 
from the present technique are compared with that from the Monte 
Carlo simulation For the range of SD of input RVs considered results 
from the simulation are almost identical to those from the present 
approximation Hence we can see that the first order perturbation 
approximation assumption for the random variables is valid foi the 
present class of pioblems. One can further conclude that neglecting of 
the higher order terms does not introduce any significant error m the 
results. 

6 5 2 Simply Supported Rectangular Plate 

Numerical results are presented for initial buckling of laminated thin 
rectangular composite plates The results are presented for a typical 
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Graphite/Epoxy system Heie we consider plates with all the four sides 
simply suppored The plates are assumed to be subjected to uniformly 
distributed compressive load form the pair of opposite X-edges 

The mean of normalized critical bucklmg load Ncr plotted against 
the aspect ratio is given in Figure 6 2 for a [ 90 °] lamina 

Figuie 6 3 shows the change in SD of normalized critical bucklmg 
load Ncr with the change m aspect ratio a/b The results are for a [ 90 °] 
lamina Plots are given for various values of SD of SD of all othei 
basic material parameters are assumed to be zero Similar plots^ for 
change in SD of i?22 / £^12 and Gn, with SD of all the other RVs kept 
zero IS given m Figures 6 4,65 and 6 6 respectively SD of Eu and 
E22 has the maximum influence on the SD of buckling load Ncr SD 
of U12 has the least mfluence on the buckling load In the case of G12 
the SD of Ncr is independent of the aspect ratio It dependent only 
on the mput SD For 1/12 except at very low aspect ratios SD of Ncr is 
independent of aspect ratio But in the case of and E22 there we 
can see abrupt jumps in the SD of Ncr as the aspect ratio of the plate 
is gradually mcreased In the first case SD of Ncr shows an mcreasing 
trend as the aspect ratio is increased. But in the case of E22 the trend is 
reversed and the SD of Ncr decreases for higher aspect ratios Further 
m both the cases for low levels of SD of input RVs, the SD of Ncr 
seems to be independent of the change in aspect ratio It may be noted 
here that En is a fiber dominated property, where as E22 is a matrix 
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dominated property 

The ratio of SD to mean is assumed to be same for all the input RVs 
The mean of noimabzed critical bucklmg load plotted against the aspect 
ratio, for laminate configurations [0°] , [90°/0°]j and [0“/90“]j are given 
in Figure 6 7, 6 8 and 6 9 respectively These indicate a trend similar 
to the results from a conventional deterministic analysis Laminate 
configurations with the outer plies having higher stiffness in the load 
direction gives higher buckling loads 

Figure 6 10 show the change of SD of Ncr with the aspect ratio for 
[90°] lamina Different plots are give to show the change in SD of Ncr 
with the change in SD of input RVs bf- The SD of Ncr has similar 
qualitative behavior with variahon in aspect ratio and m, as the mean 
of Ncr Increased randomness in the input is reflected in the increased 
level of SD of Ncr for the plate Figure 6 11 show the change of SD 
of Ncr with the aspect ratio for [90°/0°]j laminate The general trend 
is similar in nature to the [90°] lammate Figures 6 12 and 6 13 show 
corresponding results for [0°] and (0°/90°]5 lammate configurations The 
trend in the plot is similar to the plot of mean of Ncr We can observe 
that the Ncr is less and less affected by the change in aspect ratio as 
the SD of input RVs is decreased At 1 percent SD of the input RVs the 
Ncr IS virtually unaffected by the change in aspect raho of the plate 
Thus, we may be able to set a limit for the SD of mput RVs below 
which their mfluence on the aspect ratio can be neglected 
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We can see that Figure 6 10 — where al] the jnput RVs are assumed 
to have the same SD — is a combination of Figures 6 3, 6 4, 6 5 and 6 6 
The nature of variation of each one of the four basic random variables 
is different The combined nature of vaiiation of the SD of Ncr will 
strongly depend on the magnitude and nature of contribution of the 
individual RVs 

The mean values of normalized critical buckling load Ncr for for 
plates with various stackling sequences, having aspect ratio 1 0, are 
given in Table 6 1 All the plates are simply supported on the four 
edges The results for laminates having stacking sequences [45°] , 
[45°/0°]j / [0“/45°/90V-45°], and [-45V45°]j in the form of plots of SD 
of normalized critical buckling load against the SD of input RVs are 
given in Figures 6 14, 6 15, 6 16 and 6 17 respectively The results are 
for square plates, when all the input RVs are assumed to have the same 
SD to mean ratio The maximum bucklmg load is given by the [45°] 
lamina. We can see that, in the case of laminate the 0° layers 

don't have much effect on the SD of Ncr. The results in the case of 
[0'’/45“/90°/-45“]^ are similar to that for [0°] lamina given in Figure 6 12 
earlier, agam indicating the overall dominatmg influence of outer layers 

Figure 6 18 gives the SD of Ncr against SD of input RVs plot for a 
[60“] square lamina 

Now we examme some results for a typical Glass/Epoxy composite 
Figure 6 19 shows the normahzed mean buckling load plotted against 
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the aspect ratio The SD of normalized buckling load plotted against 
the aspect raho is given in Figure 6 20 Both the results are for a [90“] 
rectangular lamina with all the edges simply supported The general 
trend foi Ncr is similar in the two material systems Howevei, the 
Graphite/Epoxy mean critical buckling load plot given by Figure 6 2 
shows lower buckling loads in comparison with that for Glass/Epoxy 
composite given by Figure 6 19 The correspondmg plots for SD of 
Ncr also shows higher values for Graphite/Epoxy system Here also, as 
observed earlier at low levels of SD of input RVs change m the aspect 
raho do not have any effect on the SD of Ncr, except at very low values 
of aspect ratio 

6 5 3 Clamped Rectangular Plate 

In this section we consider thin rectangular plates clamped along all the 
four edges, subjected to a uniformly distributed mplane buckling load 
along the pair of opposite X-edges 

Table 6 2 gives the mean values of Ncr for square plates clamped on 
all the four edges having various stacking sequences Results for SD of 
Ncr for laminates having stacking sequences [90°] , [0°] , [90“/0“]jj and 
[0"/90°]„ are given in Figures 6 21, 6 22, 6.23 and 6 24 respectively Here 
the SD of Ncr is ploted against the SD of input RVs All the input 
RVs are assumed to have the same SD to mean ratio We can clearly 
see that those configurations with outer laminas having higher bending 
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stiffness give higher values of SD of Ncr 

Figures 6 25 , 6 26 , 6 27 , 6 28 and 6 29 give SD of Ncr against SD 
of input RV plots for [45“] , (4570“]^ , [-457457 , [07457907-457 
and [60°] laminate configurations In all these cases we can see the 
dominating influence of the properties of the outer layer on the SD of 
Ncr 


6 6 SUMMARY 

The peituibation technique described in Chapter 3 is applied to the 
problem of buckling of thin rectangular composite plates with random- 
ness in material propeities The validation of the technique for the 
present set of problems is done by comparing the results obtained from 
an independent Monte Cailo simulation with the same boundary con- 
ditions Stability chaiacteiistics for various laminate configurations with 
change in slacklmg sequence, aspect ratio and boundary conditions are 
presented 

The analysis shows that the fust order perturbation approximation 
foi RVs IS valid for the present class ot problems In practice the nature 
of variation of SD of critical buckling load Ncr will depend on the the 
magnitude and the nature of contribution of the individual RVs The 
results show that below a limit of SD (around 1 percent of mean) of 
input RVs the aspect ratio doesn't have much influence on the buckling 
chaiacteristics Fuithei it is observed that SD of critical buckling load 
IS strongly dependent on the buckling mode shape, if we consider 
variation of RVs -Eu and E 22 alone This indicates that for two different 
buckling modes having similar mean critical buckling load values, the 
SD of buckling load can have widely differing values The outer plies in 
laminates, in general have a dominabng influence on the SD of critical 
buckling load 
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Figure 6 1 [90°] laminate, Normalized SD of buckling load Ncr, plotted against 
the aspect ratio For various input parameter SDs Compared with results from 
a Monte Carlo simulation All sides simply supported Graphite/Epoxy 



Figure 62 Normalized mean buckling load Ncr, plotted against the aspect 
ratio a/b [90°] laminate All sides simply supported Graphite/ Epoxy 
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Figure 6 5 [90°] laminate, all sides simply supported Sensitivity of SD of 

normalized buckling load Ncr to SD of input RV x /12 SD of all other input RVs 
kept zero Graphite/Epoxy 



Figure 6 6. [90“] laminate, all sides simply supported Sensitivity of SD of 

normalized buckling load Ncr to SD of input RV G 12 SD of all other input RVs 
kept zero Graphtte/Epoxy 
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Figure 6.7. Normalized mean buckling load Ncr, plotted against the AR [0°] 
laminate All sides simply supported Graphite/ Epoxy 



Figure 6.8 Normalized mean buckling load Ncr, plotted against the AR 
[90°/0°], laminate All sides simply supported Graphite/Epoxy 
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Figure 613 SD of normalized buckling load Ncr, plotted against the aspect 
ratio a/b [0°/90°]j laminate For various input parameter SDs All sides simply 
supported Graph ite/Epoxy 
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Figure 6.14 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [45“] laminate, AR=1 All sides simply supported GraphIte/Epoxy 
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Figure 6 15 SD of normalized buckling load Ncr, plotted against the 
SD of input RVs [45°/0°J^ laminate, AR=1 All sides simply supported 
Graphite/Epoxy 
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Figure 6 16' SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [0°/45"/90°/— 45°]3 laminate, AR=1 All sides simply supported 
Graphite/ Epoxy 
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Figure 6.17 SD of normalized buckling load Ncr, plotted against the SD 
of input RVs [~45°/45°], laminate, AR=1 AH sides simply supported 
Graphite/ Epoxy 
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Figure 6 18' SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [60°]^ laminate, AR=1 All sides simply supported Graphite/Epoxy 
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Figure 619 Test case Glass/Epoxy composite Normalized mean buckling 
load Ncr, plotted against the aspect ratio a/b [90®] laminate All sides simply 
supported 
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Figure 6 20: Test case Glass/Epoxy composite [90°] laminate, Normalized 
SD of buckling load Ncr, plotted against the aspect ratio. For various input 
parameter SDs All sides simply supported 
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Figure 6 21 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [90°] laminate, AR=1 All sides clamped Graphite/Epoxy 
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Figure 6 22' SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [0°] laminate, AR=1 All sides clamped Graphite/ Epoxy 
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Figure 6.23 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [90“/0°]^ laminate, AR=1 All sides clamped Graphite/Epoxy 



Figure 6.24. SD of normalized buckling load Ncr, plotted against the SD of 
input RVs (0“/90°], laminate, AR=1 All sides clamped Graphite/Epoxy 
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Figure 6 25 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [45“] laminate. AR=1 All sides clamped Graphite/Epoxy 
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Figure 6 26 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [45°/0°], laminate, AR=1 All sides clamped Graph ite/Epoxy 
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Figure 6 27 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [— 45°/45“]j laminate, AR=1 All sides clamped Graph ite/Epoxy 
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Figure 6 28 SD of normalized buckling load Ncr, plotted against the SD 
of input RVs [0'’/45°/90°/— 45°]j laminate, AR=1 All sides clamped 
Graphite/Epoxy 
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Figure 6 29 SD of normalized buckling load Ncr, plotted against the SD of 
input RVs [60°] laminate, AR=1 All sides clamped Graphite/Epoxy 


Chapter 7 


COUPLING OF BASIC MATERIAL PROPERTIES 

7 1 INTRODUCTION 

Randomness m composites can be the result of various factors like 
variation m fiber orientation, individual lamina thickness, fiber volume 
fraction, void fraction etc Analysis of composites is normally based on 
the basic material properties like elastic modulus, Poisson's ratio etc. In 
this chapter we try to accommodate the coupling between basic material 
properties in the analysis by considering the fiber volume fraction as 
the primary random variable 

7 2 ANALYSIS 

One of the basic simplifying assumptions used in the perturbahon 
formulation presented in Chapter 3 is the independence of basic random 
variables. The random variables on which the formulation is based are 
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Gx2 


GfGm 

GjVm 4 * GjnVf 


(7 5) 


Assuming that the void content in the composite material is nil, the 
matrix volume fraction can be expressed in terms of fiber volume frachon 
[41] as, = (1 — V/) The above expressions for basic material 

parameters can be rewritten as 


•^11 = ~ Em) 


(7 6) 


E22 


EfEm 

Ef~VjiEf-Em) 


(7 7) 


1^12 = + -l^m) 


(7 8) 


G\2 


GfGm 

Gf ~ Vf{Gf ~ Gm) 


(7 9) 


Now from Chapter 3 we have the expressions for the reduced stiffness 
matrix terms Q,_, m the infinite series form as 


<312 = 

Qee ~ Gi 2 

Q44 = ^^23 


Q55 = G\2 


(710) 
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If we substitute Equations 7 6, 7 7, 7 8 and 7 9 in Equations 7 10 and 
write m a compact form, the following expression is obtained 

Ov = (711) 

Here are constants and are known functions of n, and c/ aie 

the basic fiber and matrix material properties, given by E/j Gf, Uf, E^, 
Gjji and 9t}{n) are known functions of n 

Consider the case when we study the effect on the system due 
to material parameter randomness If the fiber volume fraction Vf' is 
considered as the basic random variable, the coupling between vari- 
ous basic material properties such as elastic modulus, Poisson's ratio, 
shear modulus etc can be taken care of Thus from Equation 7 11 the 
expectation of Qj; can be found as. 

(712) 

Similarly expressions for other related quantities can also be developed 
Hence the expectation of [A^j], [B^j] and [J9,_,] matrices can be expressed 
again in the form 

= Er^®[(5„)J(/>i-Aw) 

m,] = 5Et'lTs[(5.,)J(M -/•;-,) 

= 5E^!T^[(^v)J('‘i-'>t.) (713) 

Similarly, derivatives of these terms as given in Chapter 3 can be 
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developed m the following form 


aEiQ„] 

86 ," 


U-1 


(714) 


Here /,_,j are constants and e,,i,i(n) are known functions of n, and Ck are 
the basic fiber and matrix material properties, given by Ej, Gf, vj, E^, 
G„i and are known functions of n 

From this, similar derivative forms can be developed for [A^j], [B^] 
and also 


dE[A,] 

dbj^ 

dBlB,,] 

dtp 

dEjn,,] 

dbf 


E JViajf 

k=l 




(kk — hk-i) 



iVjfay 

K=l 


iV/ay 

A:=l 


dbp 

d6f 


(hl-hU) 


(hl-hU) 


(715) 


Thus with the help of the above expressions the stahstics of the un- 
known response can be found out, following the procedure outlined 
m Chapter 3, provided we know the the statistics of the fiber volume 
fraction 


7 3 SUMMARY 

A method for taking into account the effect of mutual couphng of the 
basic material parameters like Bu, E 22 / 1^12 and G 12 m the analysis of 
composites is presented This method can be used to solve problems 
similar to those discussed in the previous chapters Thus the effect of 
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coupling of basic material parameters can be taken care of through the 
fiber volume fraction V) 



Chapter 8 


GENERAL OBSERVATIONS AND CONCLUSIONS 

8 1 INTRODUCTION 

In chapters 1 to 6, we have discussed the importance of the present prob- 
lem, the related work already done in this field and general formulation 
of the problem Specific examples with results for problems related to 
static deflection, free vibration and initial buckling have been discussed 
In this chapter we summarize some of the general observations and 
conclusions 

8 2 CONCLUSIONS 

From the results discussed m m the previous chapters it can be seen that 
the effect of randomness on the response of composites is a complex 
phenomenon dependent on a large number of parameters The general 
conclusions that can be derived from the present study are 
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• The first order perturbation approxmiation gives sufficiently ac- 
curate results for the present class of problems, for the range of 
variation m the input SD considered 

• At low levels of SD of input material parameters, change in fiber 
orientahon, aspect ratio etc do not have much influence on the 
output SDs Higher levels of SD of mput RVs show progressively 
higher influence on the response 

• The nature and magnitude of influence of each one of the primary 
random variables on the response SD is different 

% In the case of unidirectional laminae, SD of En has the maximum 
influence on the SD of the response 

• Suffer laminate configurations show lesser SD of deflection and 
higher SD of natural frequencies — if other parameters remain 
unchanged 

• The mode shape, in the case of natuial frequencies and buckling 
loads has a dominating influence on the nature and magnitude 
of variation of the response SD with the SD of input material 
parameters 

• In general if we consider bending mode, the outer plies of plate 
have an overall dominating mfluence on the SD of response. 
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8 3 LIMITATIONS 

The mam limitations of the present study can be summed up as follows 

• The technique cannot be applied to problems where the random- 
ness of the variables involved are very large, compared with the 
mean values 

• The technique is difficult to apply in its present form, to problems 
with complicated shapes and/or boundary conditions 

• The study is confined to the case of thin plates with small defor- 
mations 

• Classical laminated plate theory is used throughout, neglechng 
transverse shear effects 

8 4 SUGGESTIONS 

The followmg problems can be attempted as an extension of the present 

work 

• A method to handle coupled systems of equations, in the case of 
general composites 

• Implementation of the present technique using finite element 
method, by which the effects due to randomness m system pa- 
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ra meters can be taken care of at the element level itselt With this 
complex structural shapes can be handled easily 

• A technique to accommodate correlated basic random variables 

• Pioblems with randomness in applied static or dynamic loading 

• Problems where geometric parameters and boundary conditions 
are uncertain 
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’PENDIX A 


.SPONISE RELATIONS 


I bending stiffness matrix terms are expressed in Chapter 3 in compact form 
] as 


A, = «3., ),('■?- A*-.) 


(A I) 


IS set of equations can be expressed in terms of the primary random variables as 

^hi - kU) 

1, j = 1, 2 and 6 (A 2 } 


■re bi = £^11, 62 = E22, hs = /112, = G12 and 6$ = G23/ as defined in Chapter 3 

Equation A 2 /,j are constants, e,jt(n) the powers of 6( are functions of n — the 
unter index for the infinite series approximation of QyS, are functions of 

e fiber orientation 0 ^ for layer of the laminate, are constants and hk is the 
tckness of layer These parameters can be expressed as 

<7 = 1,2, .,7, i = 1,2 and 6 j = 1,2 and 6 (A3) 
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^111 — 1 } Cii2 — Cos { 20 )^ 

::^221 = 1 > C 222 = - Cos { 20 ), 

^qtjy ,7, z = l,2 and 6 j = 1, 2 and 6 (A 4) 

n — 3 ^ _ d 

^111 ■" gi ^122 gi 

0211 = 2 ’ *^222 = “ 2 > 

f,j, 1 = 1,2 and 6 j = 1,2 and 6 (A 5) 

/n = 1 ) /i2 = 1 ) 

/21 = 1 ) /22 = 1 ) 


Cy/, 1 = 1,2, ,5, 1 = 1,2 and 6 j = 1,2 and 6 (A 6 ) 

eiu = (1 ~ ^)) 6112 = 2 n, 

6221 = ~n, e 222 = (n + l), 


Specific terms m Dtj matrix can be expressed as 

= ^Ef:: {eU(^u,). {eLiE:=i-,u {E“o/n{n(M'^”'"^’ 

^ 12 = {eU(^i 2 ,). {eLiE:_-i-,i 2 {E:=oA 2 {n(M‘^-^"^’ 





^ 2 e = ^E.^: {e;i(^ 2 .,), (ELiE:=i-. 2 e {e“ 

= ^Er;: {e:=i(^oo,). {ELiE:=i-.oa {EL/so{n(M'“'^"' 




(hi - hU) 
- hU) 
(hl-hU) 
(hi - hU) 
{hi - hU 
{hi -hi,) 


(AT) 
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A 1 Static Deflection 


All Specially Orthotropm Case 

For the simply supported case the mean deflection can be found by solving foi the 
expectation of Lmn given by Equation 4 12 

E[L„„] = + 2(£;[A2]A;;, + 2E;[A6])A^A2 + £;p22]A^ (A 8) 

Where Am = ^ and A„ = 

From Equation 4 11, knowing the transverse load coefficient terms we can find 
out the statistical chaiacteristics of a^n and hence mean of deflection from Equation 
4 8 To solve for the SD of deflection, the expectation of partial derivatives of L,n„ 
are found using Equation 4 13 

dElLmn] _ dElBn]u , . ,2 ,2 , d£;[Z?22] u 

dbf ~ dbf db^ dbj^ \ 

(A 9) 


From the expectation of derivatives defined in Equation A 9 and using Equahons 3 15 
and 3 31, we can find out variance and hence the SD of deflection The expressions 
for Dij, given by Equations A 7 can be substituted in Equations A 8 and A 9 to arrive 
at the relations that are to be solved for finding the response statishcs for the plate 
simply supported along all edges. 

For the clamped boundary condition, the final response relations that are to be 
solved are given by M x A algebraic equations given by Equation 4 20 as 



d^X.d^Xm , 

—dx 

dx^ dx^ 
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+D 


X,X„.x /; . 4D..|" f !f..x [ 

rb pii 

= / / P{^yy)^,nyndxdy 

•Jo Jo 

m = l,2, ,M and n=l,2, ,iV (A 10) 


The expectation and expectation of derivatives of a,j can be obtained from Equation 
A 10 Using these results mean and SD of deflection can be found out from equabons 
3 15, 3 31 and 4 18 Expressions to be solved can be found by substituting Equations 
A 7 in A 10 


A 1 2 Midplane Symmetric Case 


Here the response relations are given by Equation 4 25 


r d?x, d^x, 


dx^ d.^ 

Jo dx^ dy2 dx^ dy^ 


"dx f YjYndy 

Jo 

dy^ 


+Di2 

.r. , f^<BY,SYn, , rdX,dX^ 

+2Dia 


■>> dYj dY„ 

I “Lin; I 

'o dx dx Jo dy dy 


dy"^ dy^ 

Jo dx^ dx Jo ^ dy Jo dx dx^ Jo dy 


./ 

ion f r V I J. r^dX„,_ 

+2^20 ^ ^dxj^ dy+j^ X, d dy 

n d 

pip^yV^^m^ndxdy 

_ 

m = l,2, ,M and n = l,2, ,N 


dy 


(All) 


This is of the form of Equation 3 2 Here loading p{x^y), Aj and all other terms 
except tty are known Thus Equation A 11 can be solved for the expectation of Oy and 
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expectation of derivatives of a,j, by using transformations of the form of Equations 
3 10, 3 11 and 3 13 SD of (i,j and hence that of the deflection w(^x,y^ can be found 
using the relation given by Equation 3 15 and 3 31 The expressions to be solved can 
be found by substituting Equations A 7 in the M x N linear simultaneous equations 
given by Equation All 


A 2 Natural Frequencies 


A 2 1 Specially Orthotropic Case 


For this case the final response relations to be solved are given by Equation 5 10 
Thus we have the following M x N homogeneous algebraic equations 




-\-D\2 


0 da;2 

d?X, 

d\/ 


i 


r ^ d^X, (^Xm , cPY, , 


f,i rb (fiy (fly r 

+D22j^ x.x^dx + 

-p^mnj^ X.Xmdx j^YjYndy^a,j = 0 
m = l,2, ,M and n = l,2, ,N 


di/ 

•'dX,dXm, r^dY,dY„ 


dx dx 


■dx 


j 

Jo 


dy dy 


dy 


(A 12) 


For a non-tnvial solution, the coefficients associated with must go to zero for 
all values of i and j The resulting eigenvalue problem can be solved for mean of 

and SD of by transformations of the form of Equations 3 10 and 3 11, along 

/ 

with Equations 3 15 and 3 31 The final exact relations can be found by subshtuting 
Equations A 7 in A 12 given above. 
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A 2 2 Midplane Symmetric Case 


The final response relations for this case contains terms containing and D 21 , m 
addition to those in Equation A 12 


j /" drX,drX,„ 
S.=i^,=i dx^ da2 




r" f\r , , r V f\, d% ^ 


X.X„dxl -^-^dy + iD^l 


+2£)i(j 

+21)20 




cPX, dX^ , fb _ dY„ 


di? dx 


-dx 


r dX, dX^ 

dx dx 

■'• dX, 


Jo dy ay 

f" ‘‘y- j , r j f ‘ ,/ ‘'n , ■ 

Jo Jo 


’ r .. ^ , /•'■ ^ dXm , f'‘ d?Y, dY„ , 

Jo dx ^‘'^/o dy dy-^ do; dy^ dy 

-pcal„ 1“ X,X„dx iy,Y„dy^a,, = 0 
m = l,2, ,M and n = l,2, ,fV 


(A 13) 


An eigenvalue problem is obtained from Equation A 13 by applying the condition 
for non-trivial solution This eigenvalue problem given by equating the coefficient 
term of to zero can be solved for mean and SD of natural frequencies by 
employing transformations given by Equahon 310 and 311, along with Equation 
3 15 and 3 31 Here the final expressions to be solved can be found by substituting 
Equation A 7 in Equation A 13 
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A 3 Critical Buckling Load 


A 3 1 Specially Orthotropic Case 


For the simply suppoiled boundary condition case the final ciitical buckling load 
expression is given by Equation 6 11 

TT^ r 

piim' + 2 {Di2 + 2Deo)m'^R^ + D22R‘'] (A 14) 

Mean of can be found by finding the expectation of Equation A 14 SD of Ncr can 
be found using Equation 3 15 along with 3 31 The final expressions can be found 
by substituting Equations A 7 in the above 

For the clamped boundary condition case the final equations to be solved are 
given by Equation 6 19 




ax 




+D 


12 


dx^ dx^ 

■ r d?x, , (PY„ , r ^ , ' 

Jo fix2 Jo dt/2 dx^ Jo dt/2 

«/ U wiX/ «/ U ^\j LLiXf t/ U iJbij 


lX,dX^ \ _n 


dx dx Jo dy dy 


j 

Jo 


'-dy 


+A^, 


■I 


771=1,2, ,M and n = l,2, ,// 


(A 15) 


For non-trivial solution of this equahon the coefficient of must go to zero This 
gives an eigenvalue problem, which can be transformed to a form similar to Equations 
3 10 and 3 11 and solved for mean of critical buckling load Ncr Now, using Equations 
3 15 and 3.31 the SD of Ncr can be found out The final expressions to be solved 
can be found by substituting Equahons A 7 m the set of M x N equations given by 
Equation A 15. 
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A. 3. 2 Mi({pl;)iu! Syniniotric Cnse 


lU'iv the final to bo .solved are given by Equation 6.22: 

. ,A/ ^ [ /- (/•' .V. (/■•' .V,„ rl> 

\ "i) ' '' i, 


I) 


11' 




•'<' ‘'i» f/.V (hr Jo dx dx Jo 


12/),,. 
12/ J’.'ii 


d\/ 

dxj dxj 


dy 


(" f''v ;, , r d^id^X,,, , dYi , ■ 

f .V, r,v,. 

uj/ dir J{) dx J{) dy^ dy 


HI 1,2,..., l\f and n ; 1, 2, . , . , N 


(A.16) 


Hero also the oij’onvaluo iM'<iblc‘n\ to bo solved, for finding the critical buckling load 
Mcr is obtained l»y I'niiating ll\o eoeffident term of a„„, in Equation A.16 to zero. The 
moan of Nvr is found by applying the transformations given by Equations 3.10 and 
3.11 to tlii.s I’ij’onvahu' problem. The SI) of Na' is obtained by employing Equations 
3.31 ,ind 3. IS. Tlie exact final exprc.'isions to be solved for finding the statistics of 
the erilical bueklinp, load for this case can be found by substituting Equations A.7 in 
Du; .set of AJ x N linear eqiialion.s A.16. 



APPENDIX B 


GOVERNING EQUATIONS OF ANISOTROPIC PLATES 

’llu' ivnu'inmi’ (Ujualions .is }-,iv('n in Ihe Ix'j'inning of Chaplei 4, 5 and 6 are taken 
(loin u'fc'n'iiu- | fOJ A j'.eiu'uil ileiiv.ilion foi these equations aie presented here 
\vj))) lilt' jiliysji.il cs/'JinialJon oi Die Dif/erent Icrms The co-oidinate system 
usetl IS }’;ivi‘ii I'V 1 ij'iiie ’t 1 Ihe chspliiceinenls in i, y and z directions are denoted 
by ti, II aiul ii< lespi'i lively ‘ihe mam assumptions made in deriving the governing 
equiiliims fin ii, r anil »• aie. 

1 'Ihe )>l.ile I'onsisl.s ol an aibilraiy numbci of oilhotiopic layers The mateiial 
aive.s ol the nKlivithtal layers need not coincide with the geometric axes of the 
plate. 

2 The pi, lie is Ihin. 

3, Defomuilions of the jilale u, v and w are small compared with the plate 
llneknesH. 

4. Inplane sliains are small compared to unity. 


175 



5, All lUMilincMr k-nn.s in tlie equations of motion except those involving products 
of stress and plate slopes are neglected. 

6, Tr.insverse shear idiains are negligible. 

7, Displacements u anti v are linear functions of 2 co-ordinate. 

H. MVansverse iH)rmal strain is negligible. 

9. liach |>ly obeys 1 Uioke's law. 

Id. 'riie |>la(e is of constant thickne.ss. 


Jl. Kolalojy inertia terms are lu’gligible. 

12. No ImhIv hu'ces act on the system. 

13. ' Ivans verse shear strains vanish on the 2 surfaces. 


ITom these, ll 

U' strains are given by: 




It pi ... 

de" 

(It/ 

„ On^ dv^ 

Cry 



(Ihif 



\ , ^ ^ Xv 

/V 

y ... 

cJxdy 


Here ii*' and v" art* the inid"’pl>^i''e tangenti«il displacements and w is the transverse 
displacement. 

The cqualioms of motion for the layer of laminate are given by: 
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JL t -I- ~ (a^ + a’‘ —a*—' 

Ox \ '■ 'Ox '^OiiJ yj/ \ Ox ^ dy 


( k . k k <^''A 


ox-’^Ju ) - ''•‘W 

llori' i>^ •!’ of till' individual layers Stress and moment resultants are 

d efi I u’d as: 

.ft 

(iV, , yV,y, A , j,) “ j ^[O.^^Oyy (T,^y)dz 

{Qx^Qy) = 

A 

(M,, A7„, A/,,,) j ^^{<r^:,(ry,(Tly)zdz 

Prom p;(|ualions H,?. and \l'^ Ihe following relations can be derived: 

ON,. , ON,y 


P ilj9 


i)x 


Ov 


ilNr, ON, 

'l)r Irji 


O^Mr . 

(hv^ 


i*w \7 //“'iu O'^'W d^w d'vj 

Hcluations H.4, H.S and H.6 together form the equations of motion. 

Now till* i'on.slilutive relations for the plate ate given by. 


N, ' 




VV.,; 
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'V , 
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*‘hi 

/h'.! 

/ 1 1(1 

Bu 

Bn 

Bn' 


1 

>ll 2 

M'l 

/1 20 

Bn 

Bn 

Bn 


c“ 

/lu! 

^laii 

>icfl 

Bn 

Bze 

Bqq 


e° 

'•■r.y 

fhl 

Bn 

B,o 

Dn 

Dn 

Dn 


X<r. 

Ih'i 

Bn 

Bn 

Dn 

D22 

D2, 


Xy 

/hfi 

Bn 

Bm 

Dn 

B 20 

Boo. 


. Xr-y . 


where the various stiffness matrix elements aie given by. 


(B.3) 


(B.4) 

(B.5) 

(B.6) 


(B.7) 


(B.8) 
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From Ikumlicms ll.-l, H.S, 11,6, slniin-displaccmenl and curvature-displacement rela- 
tions ,;iven by i;.|U,ilions It.l and coustitutive relations given by Equations B.7, we 
catt arrive at equations „1 itioiinn in terms of displacements for the case when any 
Iransvi'i'sc lo.ul />(•'■. ,</) .M lini', on tho plnlc ns; 


.1 




.■/•ff" 


1 1 






, (.do. t 

OXOn ()y^ Ox^ 

= 0 


(B.9) 
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-n ae, //,, 
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lir, I iAv\ 
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rV» 

O-V 
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ihH)y - (V 

fl'if) 

■ /y-., 

ihj 


(B.IO) 


,, <y'w /)'((> .,,,, . dhu d'^vj 

,, iV'rn iFu“ „„ ()V , a’u“ 


(/■')/' (‘/'o'* . (I'H/’ 



‘9'"" o ^ 


(B.ll) 


It is obscrvi'tl Dial fur a I'laurral Inminalc the governing equation in u, v and w 
are coupled, ’rhe.se e<|uations .Hinrplify for specific configuration of the stacking of 
the lamina, wlierein several terms dropout depending on the coefficient matrix of 
Equation B,7. 

In the case of free vibration analysis the general governing equations are obtained 
in the form given below. Here the load coming on the plate due to free vibration is 
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taken into aa’ounl by Ihr terms containing effective density of the lamina />. These 
terms acci>unt for the force due to inplane and transverse acceleration. 


d-a" „ IPw 

t (.'l|v t . I t- — /iii7r~r 

<hi-(hf (hi- " dx'^ 

'tn nt I ‘)f3 \ n 

Ml. fiM'.! I- “0(Ki hc-rpT - B-id-r-r 

Oxch/- dy^ 


’ 8 tr~ 


(B.12) 


.'luc I ( .d,, t I- 'I,,,--- -I- 71,107: 
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d,r;p:r - ,/h.r7;p " {Bn + 2Bm)i^ 
< hr till Oy da:'* ch'rdy 


(B.13) 


(dill (dn 

lie’ * **’('f.r'<h/ 

(die .. fdii" 


/hi 7 I ■! /h, T 2{Dn -I- -I- 


d'*w 


0\o 

' Ox“dy^ 


dr‘ ‘"d:i;**dy dxdif 


) /).... 11 /hr'. 

Of/' ( 

"..If 

(di)" (l*'ii" d“iH „ 

- "-'Mi;/" ''P “ " 


(B.14) 


Here also, dejUMuJinji', on the iieluiil laminate configuration and the terms present 
in the [/I], |/f|, [/>) matrix of lh|ualion 13.7, several terms dropout. Generally the 
inplane .sliffnes.s i.s lug, her tlum the bending stiffness. Therefore if one is interested 
in lh(? finsl few natural frt'tjuencics of the plate, the transverse oscillation frequencies 
are irnporlanl n.s coumparud to inplane oscillations. Hence one can simplify the 
equations further by ignoring the inplane acceleration teims. 

In the case of stability analysis the effect of inplane loads can be taken 
consideration by using post buckling equilibrium displacement relations, along 
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tiu' con.stilulivr ivl.ili.ms }>,ivi*n by l-quations B7. The effect due to the inplane loads 
is oxpri'ssi'd in the fnnu nf an effective lateral load. Thus the governing equations 
are obtained .ee 
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(B.16) 
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d' III 
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dj/‘ 


d'a" 

fir' 


''dt/'* 


(B.17) 


' ''‘d.r,v«r 

Modificationii and sinii*lifu'ation of the governing equations given above are often 
pottsible by liiking, into account the actual laminate configurations used. 



APPLNDIX C 


ENERGY Ml.TI lODS AND RAYLEIGH-RITZ TECHNIQUE 

C. I I lu'if'.y Mi'tluxls 

SliMiii PHtMp.v p' lit ii ‘.liiuluu’, vvhrn it is subjected to defoimations due to 

fou us .H tiiij’, (HI It I liu lol.il stiaiii uiu'ij'j slotud in a body is given by the pioduct of 
foiuu tind 1 ot n‘S| KHidiiif, displiU (.'inuiit vuitois When external foices act on a body, 
tonusj'oiuliii)', lusistivu louus aiu developed internally If the system is conservative 
the iMleiiuil iuMs|iV(‘ loites dis.ippetU and the body leluins to its unstrained state, 
upon leinoval ol the extcinal (ones. If vve considei the law of conservation of 
eiu'igy, the sum hilal ol the ivoik done by the internal foices should be equal and 
op)U).sil{- to the wtuk done by the external loads. Thus foi a system with an applied 
force of }\ muh'rp.oiiij', a viitiial displacement 61 , along the direction of the applied 
foice, the lesullui}; viilual woik done Aik, is given by 

(Cl) 
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Alsd if ('i' i i'^ flu* 1 1 


unc,.- ui cnwgy due to the vultml displacement we get 


,11' 


(MI', 


(C2) 


Woik i‘a tillin’ hv loiii’'. iU'voIojh’lI inU’iniilly 
foj(i"v ^vi‘ h.ivi' 


If IS the woik done by internal 


lIiMui' \vi' I .111 ivnli' 


-- 0 


(C 3) 


f>((l I \;)--() 

If llu’ Uil.il |»olfiiluil ,1,, {! I 1 ', horn llu’ above we can wiile’ 

(Mr,, 0 


01 


(C4) 


{C5) 


U \ r, 7r,, ■ ( oHfitani (C 6) 

II llu’ .ilnivo lomlilioii'. aic .ijijilii’tl to total potential cxpiession foi a stiuctural system, 
llu' p,nvt'inmj' difli’icnli.il ciiu.iljon!) and boundaiy conditioms can be obtained 


C2 Rayleiph Rit/ lochniquc 

Raylcip,h Rit/ li’dniiijiie makci. use of the principle of minimum potential energy 
given by luiimlion (Ml. 'nie principle slatob that of all the possible eqmllibrium states 
that satisfy bmindary conditions of a system, that which makes the total potential 
minimum lend to m.ike the system stable. 
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lo solvt' l(M llu' H'spnnsi’ ol sliiuliiu' an nppioximaling function satisfying at 
least the I’.i'oini'liu lunuulaiv miKlitions of Iho system, with unknown constants is 
assiinu'il 1 Ills t an I'c ot the loi m 

iV 

(C7) 

11 I 

wlu'ie i lUe unknown loiistanls and aie functions that satisfy at least the 
geonu’liK Ih'iuuI.uv loiuhtions ol the system This appioximating function is sub- 
stituted into tin oneie.v londilum specified by the piinciple of minimum potential 
eiu'ij’v j’lven hv llu 1 ipi.iliou f 1 'Ihe lesulling total potential eneigy expression is 
niiniiiu.’ed willi u">peil lo the lonsl.uils (/„ The niinuni7cihon piocess gives a set 
ol suuull.nu'ous ah'.ehiau equations This set of equations can be solved to obtain 
speiilu \ahu". ol i i on c pondini* lo the solution of the system equations 

When ex.ul >.olulion ot the dilfeieiilial equation of a system is difficult, Rayleigh- 
Kit/, .ij'pio.uli i.m In- usinl, lo Inui an appioxunnte solution The accuiacy of the 
solution depends on Ihe iiumbei ol leiiiis in the appioximating function and the 
luiu lions /„ ( ho'.eii, lo satisfy llu' boiiiuiaiy conditions. 
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